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Date du d^pot international (jour/mois/ann4e) 
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Date de priority (jour/mois/ann^e) 
01/07/1999 


D^posant 

COMMISSARIAT A L'ENERGIE ATGMIQUE et al. 



1 . 1! est notlfie au deposant que radnninistratibn chargee de I'examen preliminaire international a etabli le rapport 
d'examen preliminaire international pour ta demande intemationale et le lui transmet ci-joint, accompagn^, le 
cas 6ch^ant, de ces annexes. 



2. Une copie du present rapport et. le cas §ch6ant, de ses annexes est transmise au Bureau International pour, 
communication k tous les offices 6lus. 



3. Si tel ou tel office elu I'exige, le Bureau international 6tablira une traduction en lahgue anglaise du rapport {k 
Texclusion des annexes de celui-cl) et la transmettra aux offices interesses. 



4. RAPPEL 

Pour aborder la phase nationale aupr^s de chaque office 6\u, le deposant doit accomplir certains actes (depot 
de traduction et paiement des taxes nationales) dans le delai de 30 mois ^ compter de la date de priority (ou 
plus tard pour ce qui concerne certains offices) (article 39.1) (voir aussi le rappel envoye par le Bureau 
international dans le formulaire PCT/IB/301). 

Losrqu'une traduction de la demande intemationale doit etre remise a un office elu, elle doit comporter la 
traduction de toute annexe du rapport d'examen preliminaire international. II appartient au deposant d'6tablir la 
traduction en question et de la remettre directement k chaque office elu int§ress6. 

Pour plus de precisions en ce qui concerne les delais applicables et les exigences des offices elus. voir le 
Volume II du Guide du deposant du PCT. 
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1 . Le present rapport d'examen pr6liminaire international, 6tabli par I'adnninistaration chargee de I'examen preliminaire 
international, est transnnis au d^posant.conform^ment h I'article 36. 

2. Ce RAPPORT comprend 5 feuilles, y compris la presente feullle de couverture. 

H II est accompagne d'ANNEXES, c'est-a-dire de feuilles de la description, des revendications ou des dessins qui ont 
ete modifiees et qui servant de base au present rapport ou de feuilles contenant des rectifications faites aupr^s de 
radministration chargee de I'examen pr§liminaire international (voir la regie 70.16 et Tinstruction 607 des Instructions 
administratives du PCT). 

Ces annexes comprennent 5 feuilles. 



3. Le present rapport contient des indications relatives aux points suivants: 

I S Base du rapport 
Priority 

Absence de formulation d'opinion quant a la nouveaut6. 1'activit6 inventive et la 
d'appiication industrieile 

Absence d'unit6 dej'invention 

Declaration motivee selon I'article 35(2) quant a la nouveaute, ractlvite inventiv 
d'appiication industrieile; citations et explications a I'appui de cette declaration 



11 


□ 


111 




IV 
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Vlll 


□ 
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Fax: +49 89 2399 - 4465 
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Formulaire PCT/IPEA/409 (feuille de couverture) (janvier 1994) 



RAPPORT D'EXAMEN 
PRELIMINAIRE INTERNATIONAL 



Demande intemationale n** PCT/FROO/0 1 804 



1. Base du rapport 

1 . En ce qui concerne les elements de la demande Internationale (les feuilles de remplacement qui ont ete remises 
a /'office recepteur en r4ponse a une invitation faite conformement a I'article 14 sont considerees dans le present 
rapport comme "initialement d^pos^es" et ne sont pas jointes en annexe au rapport puisqu'elles ne contiennent 
pas de modifications (regies 70. 16 et 70. 17)): 

Description, pages: 

1-8,10-34 version initiale 



9.9bis,35 



re5ue(s) le 



1 6/05/2001 avec la lettre du 



11/05/2001 



Revendications, N"": 

1 -7 regueCs) le 



1 6/05/2001 avec la lettre du 



11/05/2001 



Dessins, feuilles: 
1/12-12/12 



version initiale 



2. En ce qui concerne la langue, tous les elements indiques ci-dessus 6taient k la disposition de Tadmlnistration ou 
iui ont ete remis dans la langue dans laquelle la demande Internationale a §t6 d^posee. sauf indication contratre 
donnee sous ce point. 

Ces §l§ments 6taient> la disposition de I'administration ou Iui ont ete remis dans la langue sulvante: , qui est : 

□ la langue d'une traduction remise aux fins de la recherche intemationale (selon la r§gle 23.1 (b)). 

□ la langue de publication de la demande Internationale (selon la rdgle 48.3(b)). 

□ la langue de la traduction remise aux fins de I'examen pr6liminaire Internationale (selon la regie 55.2 ou 
55.3). 

3. En ce qui concerne les sequences de nucleotides ou d'acide amines divulguees dans la demande 
Internationale (le cas echeant), I'examen preliminaire Internationale a ete effectue sur la base du listage des 
sequences : 

□ contenu dans la demande Internationale, sous forme ecrite. 

□ depose avec la demande Internationale, sous forme dechiffrable par ordinateur. 

□ remis ulterieurement a Tadministration, sous forme ecrite, 

□ remis ulterieurement a Tadministration, sous forme dechiffrable par ordinateur. 

□ La declaration, selon laquelle le listage des sequences par 6crit et fourni ulterieurement ne va pas au-del& . 
de la divulgation faite dans la demande telle que deposee, a ete fournie. 

□ La declaration, selon laquelle les informations enregistrees sous dechiffrable par ordinateur sont identiques k 
celles du listages des sequences Presente par ecrit, a ete fournie. 
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Demande intematlonale n** PCT/FROO/01 804 



4. Les modifications ont entrain^ I'annulation : 

□ de la description, pages : 

□ des revendications. n*** : 

□ des dessins, feuilles : 

5. □ Le present rapport a 6t6 formule abstraction faite (de certaines) des modifications, qui ont 6t6 consider6es 

comme allant au-dela de ('expose de t'lnventlon tei qu'il a 6X6 d^pos§, comme il est indiqu^ ci-apr^s (rdgle 
70.2(c)): 

fToute feuille de remplacement comportant des modifications de cette nature doit Stre indiqu4e au point 1 et 
annex^e au present rapport) 

6. Observations complementaires, le cas ech6ant : 



V. Declaration motivee selon I'article 35(2) quant a la nouveaute, I'activite inventive et la possibilite 
d'appiication industrielle; citations et explications a I'appui de cette declaration 

1. Declaration 

Nouveaut6 Oui : Revendications 1 -7 

Non : Revendications 

Activite inventive Oui : Revendications 1 -7 

Non : Revendications 

. Possibilite d'application industrielle Oui : Revendications 1-7. 

Non : Revendications ' - 



2. Citations et explications 
voir feuille separee 



Formulaire PCT/IPEA/409 (cadres l-VIII. feuille 2) GulHet 1998) 



RAPPORT D'EXAMEN Demande internationale PCT/FROO/01 804 

PRELIMINAIRE INTERNATIONAL - FEUILLE SEPAREE 



1 . Ilest fait reference aux documents suivants: 

Chu et al: 'First break refraction event picking using fuzzy logic systems' 
IEEE Transactions on Fuzzy Systems, US, IEEE Inc. New York, vol. 2, no. 4, 
1 novembre 1 994, images 255-266 
EP-A-0 468 229 (HNC Inc) 29 Janvier 1 992 

2. Concernant ie point V: Declaration motivee selon Tarticle 35(2) quant a la 
nouveaute, I'activite inventive et la possibilite d'applicatlon industrielle; 
citations et explications a I'appul de cette declaration 

La presente demande remplit les conditions enoncees dans V article 33 CBE, 
I'objet de la revendication 1 etant nouveau et impliquant une activite inventive, les 
raisons sont les suivantes: 

Quant a la revendication 1 : ' 
D1 decrit: 

- Systeme d'intelligence artificielle de classification d'evenements, d'objets ou de 
situations a partir de signaux et de parametres discriminants issus de modeles, 
caracterise en ce qu'il comprend (voir abrege, lignes 7-10 et page 257, colonne 
droite, lignes 8-11) ' 

- au moins une branche de traitement comprenant un systeme expert flou (voir 
page 258, colonne gauche, fig. 3) preriant une decision a partir de proprietes de 
haut niveau et de parametres discriminants de plus bas niveau extraits des 
signaux par des procedures de type traitement de signal (voir page 258, colonne 
droite, lignes 17-27 et page 265, colonne gauche, lignes 30-37), 

D1 ne decrit pas que Ie systeme expert flou est capable d'expliquer sa decision a 
I'utilisateur par I'intermediaire des regies triees par ordre d'applicabilite. 

D1 decrit seulement que Ie systeme expert flou utilise des regies pour simplifier Ie 
proces decisionnel (voir page 265, colonne gauche, lignes 43-45). 



D1: 
D2: 
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D2 (voir fig. 1 et colonne 14, lignes 35-44) montre un systeme reseau neuronal 
avec capacites expert pour trouver les donnees d'entree les plus pertinents. La 
contribution de chaque variable d'entree a la decision par une methode iterative 
est mesuree. 

Pour cette raison, D2 ne decrit pas non plus que le systeme expert flou est 
capable d'expliquer sa decision a I'utilisateur par I'intermediaire des regies triees 
par ordre d'applicabilite. 

Par consequent, I'objet de la revendication 1 implique une activite inventive et 
remplit done les conditions enoncees dans V article 33(3) CBE. 
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DESCPAMD 



Le document, reference [11] deer it un precede 
de pointage d ' arrivee premiere, qui est un precede de 
reconnaissance de . forme/ en traitement de donnees 
sismiques. Ce precede utilise ici des systemes de 
5 logique floue. Il prend en compte des attributs 
sismiques coirane caracteristiques . 

Le document reference [12] decrit un precede 
de mise en cfiuvre d'un systeme expert base sur un reseau 
neuronal a action directe, a etat continu entraine. "^Ce 
10 reseau comprend une ceuche d' elements de traitement 
d' entree' representatif s des variables d* entree, une 
ceuche d' elements de traitement de sortie 
representatif s des variables de sortie, et une eu 
plusieurs couches, d* elements de traitement caches. 

15. 

L' invention a pour objectif de pallier les 
inconvenients des systemes de 1 ■ art connu en proposant 
un nouveau systeme d' intelligence artificielle pour la 
classification, d ' evenements , d'objets eu de situations 
2 0 a partir de signaux et de paraonetres discriminants 
issus de medeles. 



EXPOSE DE INVENTION 

La presente invention concerne un syst&me 

25 d' intelligence artificielle pour la classification 
d' eveneraents , d'objets eu de situations a partir de 
signaux et de parametres discriminants issus de 
medeles, caracterise en ce qu'il comprend au moins .une 
branche de traitement comprenant un systeme expert flou 

30 (expert SEF) prenant une decision a partir de 
proprietes de haut niveau et de parametres 
discriminants de'^plus bas niveau extraits des signaux 
par des procedures de type . traitement de signal, et 
capable d'expliquer, sa decision a 1 ' utilisateur par 
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1 ' intermediaire . de regies triees par ordre 
d' applicabilite, 

Dans ce systeme expert flou on* effectue une 
descente de gradient sur les parametres : 
5 • X = y/a 

. s = ln(2a^) 

- r = ln(p) 
, . d 

avec : 

10 ■ y • position des ensembles flous des 

premisses 

• a : largeur des ensembles flous des 
. . premisses 

- p : poids des regies 

15 - d : degre d' activation de chaque classe 

pour chaque regie 
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[8] S. Muller, P. Garda, J.-D. Muller, Y- Cansi 
Seismic events discrimination by neuro- fuzzy 
merging of signal and catalogue features'' (Physics 
Chemistry of The Earth (A), vol. 24, 3, 
5 pp. 201-206, 1999) . 

[9] B. T. W. Cheng, D. B. Goldgof, L. O. Hall, *'Fast' 
fuzzy clustering" (Fuzzy Sets and Systems 93, 49- 
56, 1998) . ' 
[10] A. Klaassen, X. Driancourt', S. Muller, J.-D. 
3_0 Muller, - ^"Classifying regional seismic signals 

using TD3SrN-alike neural networks" (International 
Conference On Artificial Neural Networks '98, 
Skovde, Suede, 2-4 septembre 1998). 
[11] Chu et . al: "First break refraction event" 'picking 
15 using . fuzzy logic systems" (IEEE Transactions on 

Fuzzy Systems, US, IEEE Inc New York , vol. 2, 
N® 4/ 1 novembre 1994, pages 255-266). 
[12] EP-A- 0 468 229. . 
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REVEITOICATIONS 

1. Systeme d' intelligence artificielle de 
classification d' evenements, d'bbjets ou de situations 

5 a par'tir de signaux et de parametres discriminants 
issus de niodeles, caract^rise en ce qu'il comprend au 
moins une branche de traitement comprenant un systeme 
expert flou (SEF) prenant une decision a partir de 
propri^tes de haut niveau et de parametres 
10 discriminants de plus has niveau extraits des signaux 
par des procedures de type traitement de signal, et 
capable d ' expliquer sa decision a 1 ' utilisateur par • 
I'interra^diaire de regies triees par ordre 
d' applicabilite . • 

15 , 

2. Systeme selon la revendication 1, dans 
lequel, dans le systeme expert flou (iSEF)' on effectue 
line descente de gradient sur les parametres : 

• X = y/a 
20 ■ s = ln(2a^) 

- r = ln(p) 
■ . d 

avec : 

y : position des ensembles flous des 
25 premisses 

a : largeur - des ensembles flous des 
■ premisses- 

p : poids des regies ^ 
. d : degre d' activation de chaque . classe 
3 0 pour chaque regie 

3- Systeme selon la revendication 1, qui est 
un systeme multi-experts constitue d'au moins deux 
branches de traitement independantes, se configurant 
35 automatiquement par apprentissage statistique sur des 
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bases de donnees, ayant des proprietes particulieres et 
fusionnees par un systeme decisiormel de haut niveau. 

4. Systeme selon la revendication 3, dans 
lequel une branche comprend un classif ieur neuro-f lou 
(CNF) prenant ses decisions a partir de proprietes de 
haut niveau et de param^tres discriminants de plus bas 
niveau extraits des signaiix par des procedures de type 
traitement de signal. 



10 



5. Systeme selon I'une quelconq[ue des 
revendications 3 ou 4, dans lequel une branche comprend 
un reseau neuronal a connexions locales et poids 
partages (TDNN) constitues de bancs de f litres 
1*5 adaptatifs non lineaires, extrayant lui-meme des • 
informations discriminantes de representation temps- 
frequence des signaux "correspondants . 

6.. Systeme selon I'une quelconque des 
20 revendications pr^cedentes, qui , est xin syst^e de 
classification d'evenements geophysiques . 

7. Systeme selon la revendication 6, dans 
lequel les proprietes de haut niveau sent la 
25 localisation, la magnitude, I'heure et la date. . 
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The analysis of the signals is carried out using a 
global system, illustrated in Figure 3. 

Document reference [11] describes a method of 
first arrival recording, which is a method of shape 
5 recognition, by seismic data processing. This method 
uses fuzzy logic systems here. It takes into account 
seismic features as characteristics. 

Document reference [12] describes a method for 
implementing an expert system based on a direct action 
10 neural network^ with driven .. continuous state. This 
network comprises a layer of input processing elements 
■ representative .of the input variables, a layer of 
output processing elements representative of output 
variables, and one or several layers of hidden 
15 processing elements.; 

The aim of the invention is to compensate for the 
inconveniences of prior art systems by proposing a new 
system of artificial intelligence for classifying 
events, objects or situations from signals and 
20 discriminant parameters produced by models. 

Description of the invention 

The present invention relates to a system of 
artificial intelligence for the classification of 

25 events, objects or situations from signals and 
discriminant parameters produced by models,- 
characterised ; in that it • comprises at least one 
processing branch comprising a fuzzy expert system (FES 
expert) taking a decision according to high level 

30 properties and lower level discriminant parameters 
extracted from signals by signal processing type 
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procedures, and capable of . explaining its decision to 
the user through the intermediary of rules selected by. 
order of applicability. 

In this fuzzy expert system a gradient decrease is 
carried out on the parameters: 

• X = y/a 

• s = ln/2a^) 

• r = In (p) 

• d 

with: 

• y: position of fuzzy sets of premises 

• cr: width of fuzzy sets of premises 

• p: weight of rules 

• d: degree .of activation of each class for each 

rule. 
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western Europe" (Tectonophysics, 179, pp. 27- 

. 53, 1990) 

[8] S.Muller, P. Garda, J-D Muller, Y. Cansi 
^^Seismic events discrimination by neuro-fuzzy merging 
5 of signal and catalogue features" (Physics Cehmistry of 
The. Earth (A), vol. 24, N''3, pp. 201-206, 1999). 

[9] B. T. W. Cheng, D. B. Goldof, L.O. Hall, „Fast 
fuzzy clustering'^ (Fuzzy Sets and Systems 93, 49-56-, 
1998). 

10 [10] A.Klaassen, X. Driancourt, S. Muller, J.-D. 

Muller, "'Classifying regional seismic singals using 
TDNN-alike neutral networks" (International Conference 
- On Artificial Neural Networks '98), Skovde, Sweden, 2-4 
September 1998) . 

15 [11] Chu et al: "First ' break refraction event 

picking using fuzzy logic systems" (IEEE Transactions 
on Fuzzy Systems, US, IEEE Inc. New York, vol. 2, No. 
4, 1 November 1994, pages 255-266). 
[12] EP-A-0 4 68 229. 
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' CLAIMS 

1. System of artificial intelligence for 
classification of events, objects or situations from 
signals and from discriminant- parameters produced by 
models, characterised in that it comprises at least one 
processing- branch comprising a fuzzy expert system 
(FES) taking a decision according to high level 
properties and discriminant parameters of lower level 
extracted from signals by signal processing type 
procedures, and capable of explaining its decision to' 
the user through the intermediary of rules selected by 
order of applicability. 

2. System according to claim 1 in- which, in the 
fuzzy expert system (FES), a gradient decrease is 
carried out on the parameters: 

• X = y/a 

• s =' ln/2a^) 
r = In (p) 

• d 

.with: • . . ■ ^ 

• y: position of fuzzy sets of premises 

• a: width of fuzzy sets of premises 

• 'p: weights of rules 

• d: degree of activation of each class for 
each rule. 

3. System according to claim 1, which is a multi* 
expert system constituted of at least two independent 
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processing branches, organising themselves 

automatically through statistical learning on data 
bases, having particular properties and merged by a 
high level decisional • system. 

4. System according to claim 3, in which one 
branch corresponds to a neuro-fuzzy classifier (NFC) 
taking its decisions based on high level properties and 
lower level discriminant parameters extracted from 
signals by signal processing type procedures. 

5. System according to one or the other of claims 
3 or 4, in which one branch comprises a neural network 
with local connections and shared weights (TDNN) 
constituted" of banks of non-linear adaptable filters, 
itself extracting .discriminant information for time- 
frequency representation of the corresponding signals . 

6. System according to any one of the preceding 
claims, which is a system for classification -of 
geophysical events. 

'7. System according to claim 6, in which the ^high 
level properties are the localisation, the magnitude, 
the time and the date. 
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3, rue Du Docteur Lancereaux 

F-75008 Paris 

FRANCE 



Reference du dossier du dSposant ou du mandataire 
B 13340.3 DB 



NOTIFICATION IMPORTANTE 



Oemande Internationale no 
PCT/FROO/01804 



Date du depot international (jour/mois/ann6e) 
28 juin 2000 (28.06.00) 



Date de publication Internationale (ipur/mois/annie) 

Pas encore publi6e 



Date de priority tiour/mots/ann6e) 

01 juillet 1999 (01.07.99) 



Deposant 

COMMISSARIAT A L'ENERGIE ATOMIQUE etc 



1 . La date de reception (sauf lorsque les lettres "NR" figurent dans la colonne de droite) par le Bureau international du ou des 
documents de prionte correspondant k la ou aux dennandes §numerees ci-apres est notifi6e au deposant Sauf indication 
contraire consistent en un asterisque figurant d c6t6 d'une date de reception, ou les lettres "NR". dans la colonne de droite. 
le document de priorite en question a 6te presente ou transmis au Bureau international d'une maniere conforme d la 
regie 17.1. a) ou b). 

2. Ce formulaire met a jour et remplace toute notification relative a la presentation ou a la transmission du document de priority 
qui a 6te envoyee prec6demment 

3. Un asterisque(*) figurant a cote d'une date de reception dans la colonne de droite signale un document de priorite pr6sent6 
ou transmis au Bureau international mais de maniere non conforme d la regie 17.1.a) ou b). Dans ce cas, Pattention du 
deposant est appetee sur la regie 17.1.c) qui stipule qu'aucun office designe ne peut decider de ne pas tenir compte de la 
revendication de priorite avant d'avoir donpe au deposant la possibilite de remettre le document de priority dans un d6lai 
raisonnable en I'espece. 

4. Les lettres "NR" figurant dans la colonne de droite signalent un document de priorite que le Bureau international n'a pas 
repu ou que le deposant n'a pas demande a I'office recepteur de preparer et de transmettre au Bureau international, 
conformement a la rfegle 17.1.a) ou b), respectivement. Dans ce cas, Tattention du deposant est appelee sur la regie 17.1.c) 
qui stipule qu'aucun office d6signe ne peut decider de ne pas tenir compte de la revendication de pnorit6 avant d'avoir donn6 
au deposant la possibilite de remettre le document de priorite dans un delai raisonnable en I'espece. 



Date de priority Demande de priorite n 

01 juil 1999 (01.07.99) 99 08472 



PaY<; office regional ou 
office rece nteur selon le PCT 

FR 



Date de reception du 
document de priority 

02 aoOt 2000 (02.08.00) 



Bureau international de TOMPI 
34, chemin des Colombettes 
1211 Gendve 20, Suisse 



no de t6lecopieur (41-22) 740.14.35 



Fonctionnaire autorisex 



Simin Bahartou 

no de t6l6phone (41-22) 338.83.38 



Formulaire PCT/IB/304 (juillet 1998) 



003495600 



TRAITE D 



OPERATION EN MATIERE^ BREVETS 

Exp§diteur: le BUREAU INTERNATIONAL 



WO 01/03064 
PCT/FROO/0180 



PCT 



AVIS INFORMANT LE DEPOSANT DE LA 
COMMUNICATION DE LA DEMANDE 
INTERNATIONALE AUX OFFICES DESIGNES 

(rdgle 47.1. c), premiere phrase, du PCT) 



Date d'expMition (jour/mois/ann^e) 
11 janvier2001 (11.01.01) 



Destinataire: 
LEHU, Jean 
Brevatome 

3, rue du Docteur Lancereaux 
F-75008 Paris 
FRANCE 



BREVATOME 



1 9 JAN. 2001 



3, rue du Doaeur Lancereaux 
75O0 8 PARIS 



R6f6rence du dossier du deposant ou du mandataire 
B 13340.3 DB 


AVIS IMPORTANT 


Demande intemationale no 
PCT/FROO/01804 


Date du d6pdt international (jour/mots/a nn^e) 
28 juin 2000 (28.06.00) 


Date de priority Uour/nnois/ann6e) 
01 juillet 1999 (01.07.99) 


D6posant ' 

COMMISSARIAT A L'ENERGIE ATOMIQUE etc 



1. II est notifid par la pr^sehte qu'a la date indiqu^e ci-dessus comma date d'exp^dition de cet avis, le Bureau international a 
communique, comma le pr6voit rarticle 20, ta demande intemationale aux offices d^signes suivants: 

US ^ 

Conformement S la rdgle 47.1. c). troisidme phrase, ces offices acceptent le present avis comme preuye dSterminante 
du fait que la communication de la demande Internationale a bien eu lieu d la date d'expedition indiquee plus haut et le 
deposant n'est pas tenu de remettre de copie de la demande intemationale d I'office ou aux offices d^signes. 

2. Les offices designes suivants ont renonce a I'exigence selon laquelle cette communication doit etre effectu(§e d cette date: 

CA,EP,JP 

La communication sera effectuSe seulement sur demande de ces offices. De plus, le deposant n'est pas tenu de remettre 
de copie de la demande Internationale aux offices en question (r^gle 49.1 )a-bis)). 

3. Le present avis est accompagn^ d'une copie de la demande intemationale publi6e par le Bureau international le 
1 1 janvier 2001 (1 1.01.01) sous le num6ro WO 01/03064 

RAPPEL CONCERNANT LE CHAPITRE II (article 31 .2)a) et regie 54.2) 

Si le deposant souhaite reporter Touverture de la phase natipnale jusqu'd 30 mois (ou plus pour ce qui concerne certains 
offices) a compter de la date de priorite, la demande d'examen pr6iiminaire intemationat doit etre pr6sent6e d .. 
I'administration comp^tente chargSe de I'examen pr^liminaire international avant Texpiration d*un d6lai de 19 mois d 
compter de la date de priorite. 

II appartient exclusivement au deposant de veiller au respect du d^lai de 19 mois. 

II est a noter que seul un deposant qui est ressortissant d'un Etat contractant du PCT ii6 par le chapitre II ou qui y a son 
domicile peut presenter une demande d'examen pr^liminaire internationaL 

RAPPEL CONCERNANT L'OUVERTURE DE LA PHASE NATIONALE (article 22 ou 39.1)) 

Si le deposant souhaite que ia demande intemationale procdde en phase nationale, il doit, dans le d§lai de 20 niois ou 
de 30 mois. ou plus pour ce qui concerne certains offices, accomplir les actes mentionh^s dans ces dispositions auprds 
de chaque office design^ ou elu. 

Pour d'autres informations importantes conoernant les d^lais et les actes d accomplir pour I'ouverture de la phase 
nationale, voir I'annexe du formuiaire PCT/IB/301 (Notification de la reception de Texemplaire original) et le volume II 
du Guide du deposant du PCT. 



Bureau international de I'OMPI 


Fonctionnaire autoris^ 


34, chemin des Colombettes 


J. Zahra 


1211 Geneve 20, Suisse 


no de telecopieur (41-22) 740.14.35 


no de telephone (41-22) 338.83.38 



Formuiaire PCT/IB/308 Guillet 1996) 



3752292 



La demande d'examen preliminaire inUmational doit etre presentee directement a I administration chargee de I'examen preltmmatre 
international qui est competente ou, sipiusieurs administrations sont competentes, a I une d entre elles, au choix du deposanL U deposant 
pent indiquer le nom complet ou le code a deux lettres de eette administration au dessus de la ligne qut suit 



IPEA/, 




PCX 



CHAPITRE II 



DEMANDE D'EXAMEN PRELIMINAIRE INTERNATIONAL 

selon rarticle 3 1 du Traiti de coop^iadon en inatiere de brevets : 
soussigne lequieit que la deinaiide inteinatioiiale spdcifide d-^;ir6s &sse I'objet 
d'un examen priliiniiiaire intemational confpnnimeiit au Traitd de cooperation en matieie de brevets et 
&it flection de tous les Etats digibles sauf indication contiaire. 



Administration chargdc dc Tcxamcn preliminaire international 


Date de r6ception.de la demande d^ examen preliminaire intemational 


Cadre n» I IDENTIFICATION DE LA DEMANDE INTERNATIONALE - 


B 13340.3 DB 


Demande intcmationalc n* 
PCT/FROO/01804 


Date du depdt intemational Qour/mois/annee) 
28 juin 2000 (28.06.2000) 


Date de priorite (la plus ancienne) 
(jour/mois/emnee) 

01 juillet 1999 (01.07.1999) 



litre de P invention 

SYSTEME DMNTELLIGENCE ARTIFICIELLE POUR LA CLASSIFICATION D'EVENEMENTS. □•OBJETS OU DE 

SITUATIONS. . ' ' 



Cadre II DEPOSANT(S) 



Nom et adresse • (Som defamille suivi du prenom: pour ime perspnne morale, designation officielle 
iNom ci aarcMc J.^^^^^ 'adresse doit comprendre le code postal et le nom dupayi.) 

COMMISSARIAT A L'ENERGIE ATOMIQUE 
31-33, rue de la Federation 
75752 PARIS 15eme 
FRANCE 



n* de telephone 

01 69 08 82 96 



n* de telecopietir 

01 69 08 82 92 



n* de teteimprimeur 



Domicile (nom de PEtat)^ : 
FR 



Nationalite (nom de Pfetat) : 

FR " 

Nom et adresse : (Horn defamille sunn du prenom: pour une personne morale, designation officieUe complete. L'adresse doit comprendre le code postal et le nom du pays) 

MULLER Jean-Denis 
5 allee de la Garenne 
78120 CLAIREFONTAINE-EN-YVELYNES 
FRANCE 



Nationalite (nom de I'Etat) ; 


Domicile (nom de PEtal) : 




FR 


FR 




Nom ct adresse : (Som de famille suivi du prenom; pour une personne morale, des 


ffiation officielle compleU. I 'adresse doit comprendre le code^ 


wstal et le ttom du pays) 



MULLER CARCELES Stephanie 

5 allee de la Garenne 

78120 CLAIREFONTAINE-EN-YVELYNES 

FRANCE 



Nationalite (nom de I'^tat) : 
FR 



Domicile (nom de PEtat) : 
FR 



I I D'autres deposants sont indiques sur unc fcuillc annexe. 



Formulaire PCT/IPEA/401 (premiere fcuille) (juillet 1998; reimpression juillet 2000) 



Voir les notes relatives au formulaire de 
demande d 'examen preliminaire intemational 



Fcuille n" 



Demandc intemationalc 

PCT/FROO/01804 



Cadre n» HI MANDATAIRE OU REPRfeSENTANT COMMUN; OU ADRESSE POUR LA CORRESPONDANCE 



U pcrsohnc indiquec ci-dcssous est [x] mandateire | \ representant commun 

ct [x] a designee a une date antcricurc; elle rcprcsentc aussi le ou les diposants pour Pcxamcn prcliminaire international. 
I — I est designee par la prcsentc; toutc designation antcricure dc mandataircs ou d*un representant commuri est dc cc fait rcvoquec. 

• est disignec par U prcsente. spccialcmcnt pour la procedure devant P administration chargce dc Pexamcn prcliminaire 
I I international, en sus du ou dcs mandataircs ou du representant commun designes antdneurement. 



Norn et adresse : (Norn it famUU suivi du prenom- pour me perspme morale, daiffutiion officielk 
complete. L adresse doit eomprendre le code postal et le nam au pays.) 



LEHU Jean 



do BREVATOME 

3, rue du Docteur Lancereaux 
75008 PARIS 
FRANCE 



n* de t^liphone 

01 53 83 94 00 



n* de telecopieiir 

01 45 63 83 33 



n* de teieiix^Thmeur 



I 1 Adrcssse Dour la correspondance : cochcr cette case lorsquc aucun mandatairc ni representant commun n'cst ou n'a etc 

□ dfslTndlt que r esM^ est utUisc pour indiquer un2 adresse specials a laqucUe la correspondance doit ctrc cnvoyec. 



Cadre n» IV BASE DE L'EXAMEN PRiUMINAiRE INTERNATIONAL 



Declaration concemant les modifications :* 

1 . Lc deposant souhaite que I'cxamen prcliminaire international commence sur la base suivante : 
I I la demandc intcmation^c telle qu'clle a etc deposee initialcment 

la description Q telle qu'clle a etc deposee initialement 

I — [ telle qu'clle a etc modificc en vcrtuderarticlc 34 

Ics rcvcndications Q tcllcs qu'cllcs ont etc dcposces initialcment 

I — I tcllcs qu'cllcs ont etc modificcs en vcrtu de Particle 19 (avcc, le cas echcant, U declaration jointc aux 

' — * Riodifications) 

I I tcllcs qu'elles ont cte modificcs en vcrtu de Taiticlc 34 

lesdcssins Q tels qu'ils ont Ate deposes initialcment 

I [ . tels qu'ils ont etc modifies en vcrtu dc Particle 34 

2. rn Lc deposant souhaite que les modifications apportces aux rcvcndications en vcrtu dc Particle 19 soicnt considcrces commc 

ccartces. 

3 rn Lc deposant souhaite que lc commencement dc Pcxamcn prcliminaire international soit differc jusqu'a Pcxpiration d'un dclai 

' ' dc 20 mois a compter dc la date dc pnoritc, a moins que i' administration chargce dc Pcxamcn prclimmairc mternational nc 

rccoivc unc copic dcs modifications cffcctuccs en vcrtu dc I'articlc 19 ou unc declaration du deposant, aux termcs de laquellc 
cdui-d ne souhaite pas cfTcctuer de modifications en vcrtu dc Particle 19 (regie 69.1.d)). (Ne pas cocher cette case 
lorsque le delai vise a I 'article 19 a expire.) 
• LorsQue aucunc case n'est cochce, Pcxamcn prcliminaire international commenccra sur la base dc la demandc intemationalc telle 
Qu'ellc a etc deposee initialement ou, si P administration chargce dc Pcxamcn prcliminaire international re?oit copic dcs 
modifications apportces aux rcvcndications cn vcrtu de Particle 19 ou dcs modifications apportces a a demandc mteriiaUonde en 
vcrtu dc I'articlc 34 avant d'avoir commence a rediger une opinion ccnte ou le rapport d examen prclimmairc international, sur la 



Langue : Pcxamcn preliminaire international sera efTectue en ^^05315 , qui est 

[T] la langue dans laquellc la demandc intemationalc a etc deposee. 

I I la langue d'une traduction remise aux fins de la recherche intemationalc. 

I J la langue dc publication dc la demandc intemationalc. 

I I la langue dc la traduction (qui va ctrc) remise aux fins de Pcxamcn prcliminaire international. 



Cadre n« V ELECTION D'ETATS 



Le deposant elit tous les ttats £ligibles (c 'esi-a-dire tous les ttats qui ont eti designes et qui sont lies par le chapitre II du PCT) 
a Pexclusion dcs 6tats ci-aprcs que le deposant souhaite ne pas 6Ure : 



FonnuUire PCT/IPEA/401 (deuxic^c fcuiUc) OuiUct 1998; rcimprcssion juilic. 2000) ^/« ^^'jllZ.ZJ^Z^!::^^ 



Feuille n"* 



Demandc intemationalc 

PCT/FROO/01804 



Cadre VI BORDEREAU 



Aux fins dc rexamen prcliminaire international, Ics ildmcnts suivants, itablis dans U 
langue indiqucc au cadre n*' IV, sont joints a la presente demande d*examen : 


Reserve a T administration chargee 
de Pexamen preliminaire international 

re^u non re^u 


1 . traduction de la demande intemationale 


feuiUes 


□ 


□ 


2. modifications selon l*article 34 


feuiUes 


□ - 




3. copic (ou, si ellc est cxigcc, traduction) des 
modifications scion T article 19 


feuiUes 


□ 


□ 


4. copie .(ou, si cite est exigee, traduction) de 
la declaration selon 1* article 19 


feuiUes 


• P 




5. lettre 


feuiUes - 


□ 


□ 


6. autres pieces (preciser) 


feuiUes 


□ 


□ 



Le ou les elemenU coches ci-apres sont aussi joints a la demande d'examcn preliminaire intemaUonal : 

1 . Q] fcuiUc de calcul des taxes 4. Q explication de I'absence d'une signature 

2. I [ poiivoir distinct signc 

3. I I copie du pouvoir general; 
^— ' numero de reference, le cas © 



5. r~n listage des sequences de nucleotides ou d'aeides 
' — ' amines sous forme dechiffrable par ordinateur 



( echeant : 



6. I I z.Mtxts iA&meaXs (preciser) \ 



Cadre n« VII SIGNATURE DU PftPOSANT, DU MANDATAIRE OU DU REPRfeSENTANT COMMUN 

A cote de chaque signature, indiauer le nom du signataire et. si cela n apparaitpas clairement a la lecture de la demande d 'examen 
preliminaire international, a quel titre I'interesse signe. 



LEHU Jean 




Reserve a T administration chargee de Texamen prefiminaire international 



1 . Date effective de reception de la D EMA NDE 

D'EXAMEN PRELIMINAIRE INTERNATIONAL: 



2. Date modificc dc reception dc la demande d'examcn 
preliminaire international, en cas de CORRECTIONS 
. apportees en vertu de la regie 60. 1 .b) : 



^ I 1 La demande d^cxamcn preliminaire international a cte re^ue PLUS DE 19 I | 

3- I I mois aprcs la date dc prioritc ct les points 4 et 5 ne sont pas apphcables. i 1 



Le d6posant a ete informe 
en consequence. 



La demande d'examen preliminaire international a etc ra?uc dans le delai de 19 mois a compter de la date de prioritc, proroge 



en vertu de la regie 80.5 



I — I Bicn que la demande d'examcn preliminaire international ait cte rei?uc plus dc 19 mois aprcs la date dc prioritc. le retard 
I—I a Tarn vce est EXCUSE en vertu dc la regie 82. 



Reserve au Bureau international 



Demande d'examen preliminaire international re9ucdc T administration 
chargee dc Pcxamen preliminaire international Ic : 



Formulairc PCT/IPEA/401 (demiere feuUlc) (juillel 1998; reimpression juillct 2000) 



Voir les notes relatives au formulaire de 
demande d 'examen preliminaire international 



PCX 



CHAPITRE II 



FEUILLE DE CALCUL DES TAXES 
Annexe de la demande d'examen prtliminaire international 



Demande 
internationaie n" 


PCT/FROO/01804 


Reference du dossier du 
deposant ou du mandataire 


B 13340.3 DB 



preUminsire intematioiul « 



Timbre & dote de I'administratian chargee de rexamen 



prelinuiiaire intematioDal 



Deposant . ^ 

COMMISSARIAT A UENERGIE ATOMIQUE - MULLER Jean-Denis - 
MULLER CARCELES Stephanie 



Cakul des taxes prescrites 

1. Taxe d*examen prcliminaire 

2. Taxc de traitcmcnt (Les deposants de certains Etats 
ont droit a une reduction de 75% de la taxe de 
traitemenL Lorsque le deposant a (outous les deposants 
ont) droit a eeite reduction, le montant devantfigurer 
sous H est egal a 25% de la taxe de traiument) 

3. Total des taxes prescrites , 

Additionncr les montants portcs dans les cadres P ct H 
et inscrire le resultat dans le cadre TOTAL 



1 533 Euros 



148 Euros 



H 



1 681 Euros 



TOTAL 



Mode de paiement 

Hautorisation de debiter un compte de 
depot aupres de T administration 
chargec de Texamcn preliminaire 
international (voir plus bas) 

I I cheque 

I I mandat postal 

I I traite bancaire 



I I especes 

[ j timbres fiscaux 

I I coupons 

I I 2\xUc (preciser): 



Autorisation concemant un compte de depot (les administrations chargeesde I examen preliminaire international nepermettent pas 

Unites I 'utilisation de ce mode de paiement) 

L^administration chargee de Texamen 
preliminaire international/ 



2804.0035 



est autorisee k dcbiter mon compte de depot du total des taxes ihdiqu6 ci-dessus. 

□ (cetie case ne peut etre cochee que si les conditions relatives aux comptes de depot etablies 
par V administration chargee de 1 'examen preliminaire international le permettent) 
est autoriscc a dcbiter mon compte de depot dc tout montant manquant - ou A le 
crediter de tout excedent - dans le paiement du t otal des tays indique ci-dessus. 

17 janvier 2001 



Numcro du compte dc depot 




Date (jour/mois/annee) 



Signature 



Formulaire PCT/IPEA/401 (Annexe) Quillet 1998; rcimprcssion juillet 2000) 



Voir les notes relatives a la feuille de calcul des taxes 



l^'S^^'^i'i^!^ PUBmEEmyEimJ^W^TIUITEDE COOPERATION 
- - ' EN MATIERE DE BREVETS (P^^^ 



(19) Organisation Mondiaie de la Propria^ 
Intellectuelie 

Bureau international 

(43) Date de la publication internationaie 
11 Janvier 2001 (11.0L2001) 





PCX 



:iiniiiiiiiHiiiiiiiil 

(10) Numero de publication Internationale 

wo 01/03064 Al 



(51) Classification interaatioaale des brevets^: G06N 5/04 

(21) Numero de la demaiide tntemationale: 

PCr/FROO/01804 

(22) Date de depot international: 28 juin 2000(28.06.2000) 

(25) Langue de depot: frangais 

(26) Langue de publication: frangais 

(30) Donn^es relatives a la priorite: 

9908472 IjuiUet 1999 (01.07.1999) FR 

(71) D^posant (pour tous les Bats designes sauf US): COM- 
MISSARIAT A L'ENERGIE ATGMIQUE [FR/FRJ; 
3 1 -33, rue de la Federatioiu F-75752 Paris I5eme (FR). 



(72) Inventeurs; et 

(75) Inventeurs/Diposants (pour US seulemeni): MULLEK, 
Jean-Denis [H^/FR]; 5. aI16e de ia Gaieime, F-7isi20 
Clairefontaine-en-Yvdyn^ (FR). MULLER CARCE- 
LES, Stephanie [FR/FR]; 5, all^e de la Garenne, F-78i20 . 
C3airefontaine-en-Yvelynes (FR).' 

(74) Mandataire: LEHU, Jean; Brevatome, 3, rae du EkKteur 
Lancoieaux. F-75008 Paris (FR). 

(81) Etats designes (national); CA, JP. US. 

(84) Etats designes (regional): brevet europeen (AT, BE, CH, 
CT, DE» DK, ES, H. FR, GB, GR, IE, IT. LU. MC, NL, PT. 
SE). 

Publite: 

— Avec rapport de recherche internationaie. 

[Suite sur la page suivante] 



v.. 



(54) Tide: ARTIFICIAL INTELLIGENCE SYSTEMS FOR CLASSIFYING EVENTS, OBJECTS AND SITUAnONS 

- ^5 (54) Titre: SYSTEME D'INTELLIGENCIE ARTIFICIELLE POUR LA CLASSfflCATION D'EVENEMENTS, D'OBJETS OU 
^= DE SITUATIONS 



decision globale par 
fusion, multi-experts 



CNF 



donn6es de 
hautni^au 



discriminants 
extraits des signaux 



SEF 



donates de 
hautfweau 



discrinwiants 
extraits des signaux 



TDNN 



pointesdes 
fihasesPetS 



representations temps- 
fr^ence des signaux 



A. ..GLOBAL DECISION BY MULTKEXPERT MERGING 

B. ..NEURO^UZZY CLASSIFIER 

C. ..FUZZY EXPERT SYSTEM 
□...HIGH-LEVEL DATA 

E...D1SCRIMINANTS RETRIEVED FROM SIGNALS 



F.. .HIGH-LEVEL DATA 

G-DISCRIMINANTS RETRIEVED FROM SIGNALS 
H... NEURAL NETWORK WITH LOCAL CONNECTIONS 
I...PLOrnNG OF PHASES P AND S 
J.„TIME-FREQUENCY REPRESENTATIONS OF SIGNALS 



(57) Abstract: The inventioD concerns an artificial intelligence system for classifying events, objects and situations from signals 
\^ and discriminating parameters derived from models comprising at least a processing branch including a fuzzy expert system taking a 
O decision on the basis of hig|i level properties and low level discriminating parameters retrieved form the signals by signal processing 

procedures, and capable of explaining its decision to the user using ndes s(»ted out in order of applicability. 



(57) Abrege: La pr&ente invention conceme un systfeme d'intelligence artificielle pour la classification d'evenements, d^objets ou 
de situations h partir de signaux et de param^ties discriminants issus de mod^es comprenant au moins une branche de traitement 
comprenant un syst^me expert flou prenant ime d^ision k paitir de propri^t^ de haut niveau et de param^tres discriminants de plus 
bas niveau extraits des signaux par des (HXK^ures de type traitement de signal, et capable d'expliquer sa d&iston k Tutilisateur par 
^ Tinterm^diaire de regies triees par ordre d* applicability. 



wo 01703064 Al 
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En ce qui concerne les codes a deux lettres et autres. aBreviar 
, . , tipns,.se. r^/f[er~ omx "iilotes.explicatiyes. relatives aux . codes^et ' 
^;Payreviati6ns: ^^ic^ au dibut de'iiu^ oriUnqve^dk t 

la Gazette du PCX. " ' . ^ ; 



wo 01/03064 



PCT/FROO/01804 



1 

SYSTEME D ' INTELLIGENCE ARTIFICIELLE POUR LA 
CLASSIFICATION D ^ EVENEMENTS . OB JETS OU DE SITUATIONS 

5 

DESCRIPTION 

. DOMAINE TECHNIQUE 

La presents invention concerne un systeme 
10 d' intelligence artificielle pour la classification 
d' evenements , d'objets ou de situations a partir de 
signaux et de parametres discriminants issus de 
modeles . 

L' invention s' applique notamment a la 
15 classification d' evenements sismiques. Une telle 
classification sera consideree a titre d'exemple non 
limitatif dans la suite de la description, 

ETAT DE LA TECHNIQUE ANTERIEURE 

20 La classification automatique d' evenements sismiques 

La classification automatique d' evenements 
sismiques est un probleme relativemenc recent, puisque 
le probleme n'a veritablement ete aborde que dans les 
annees 1980. Ces travaux sont ma j oritairement orientes 

25 vers la recherche de parametres discriminants (c'est-a~ 
dire susceptibles de permettre la classification) dans 
les signaux sismiques. De nombreuses caracteristiques 
potentielles sont proposees en vue d'une future 
classification automatique. A partir de 1990, on 

3 0 commence a trouver dans les articles publies des 
tentatives de classification automatique, que ce soit 
par des techniques neuronales ou par des systemes a 
base de regies. Ces travaux cherchent a separer les 
seismes d'origine naturelle des explosions. Aucun de 
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ces articles ne traite de la discrimination 
d'effondrements miniers coups de terrain "). 

Du fait de la complexite du probleme, ces 
articles mettent clairement en evidence la necessite de 
5 mettre au point des systemes automatiques capables 
d' apprentissage . Les methodes neuronales ont done 
souvent ete proposees pour realiser la discrimination 
automatique d'evenements sismiques, mais avec des 
limitations qui sont analysees plus loin. Les modeles 
10 les plus couraimnent proposes sont les percept rons 
multi-couches, a connexions completes entre couches 
successives . 

Tous ces articles cherchent a determiner 
I'origine du seisme a partir de caracteristiques 
15 extraites uniquement des signaux. Les donnees de plus 
haut niveau (date, heure, latitude, longitude, 
magnitude...) ne sont jamais employees pour la 
classification. Or, les sismologues connaissent la 
.difficulte de discriminer les signaux sismiques 
2 0 uniquement par des traitements de bas niveau. 

Les travaux de Baumgardt et de ses 
collaborateurs decrits dans le document reference [1] 
sont sans doute ceux qui ont permis les plus grandes 
avancees dans la recherche de parametres 
25 discriminants. 

Les variations du cepstre, le cepstre d'un 
signal x etant la transformee de Fourier inverse du 
logarithme de la transformee de Fourier de x, sont 
souvent utilisees. On peut ainsi montrer que le cepstre 
30 permet de visualiser le phenomene de micro-retards 
presents dans les signaux de tirs caracterises par une 
plus large variance. Le document reference [2] note 
aussi cette propriete en , signalant neanmoins que 
1' absence de cette caracteristique ne permet aucune 
35 deduction quant a la classe de I'evenement. 
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Les rapports des amplitudes des differents 
types d'ondes peuvent aussi servir de discriminants. Le 
document reference [3] etudie toute une serie de 
rapports d' amplitudes {Pn/Lg, Pg/Lg, Lg/Rg) . Ces rapports 
5 sont annonces comme pouvant permettre une bonne 
discrimination . 

Les memes auteurs introduisent aussi les 
rapports des densites spectrales de puissance des 
differents types d'ondes detectes. Comme pour les 

10 rapports d' amplitudes, ces discriminants sont utilises 
par toutes les etudes qui recherchent des discriminants 
dans les signaux sismiques. Pour caracteriser les 
explosions, on a egalement utilise les rapports des 
densites spectrales de puissance d'un type d'onde, ici 

15 S, sur des bandes de frequences differentes, c'est-a- 
dire le rapport de la densite spectrale de puissance de 
S dans la gamme 1-2 Hz a la densite spectrale de 
puissance de cette meme phase dans la bande 7-20 Hz. Le 
rapport entre des densites spectrales de puissance de 

20 I'onde S avant et apres 10 Hz est egalement presente 
comme un bon separateur des explosions des seismes. 

Le document reference [4] remarque que le 
temps de propagation des signaux provenant d'une mine a 
un temps tsg-tpg constant pour une station 

25 d' enregistrement donnee. Ce temps de propagation est 
presente comme une caracteristique potentielle pour une 
mine, restant cependant moins fiable que les 
precedentes caracteristiques . 

Le document reference [5] propose d'utiliser 

30 la presence de I'onde de surface des seismes pour les 
discriminer des explosions nucleaires a des distances 
regionales. La caracterisation de la presence d'une 
onde de surface est effectuee indirectement en 
comparant les magnitudes mb et Mg. Pour deux evenements 

35 sismiques de meme magnitude mb, la magnitude de I'onde 
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de surface Ms est en general, plus elevee dans le cas 
d'un seisme a cause de la presence de I'onde de surface 
que dans le cas d'une explosion. En effet, cette onde 
de Rayleigh de croute intervient dans le calcul de la 
5 magnitude Mg et sa presence est subordonnee au 
phenomene de cisaillement absent dans le cas 
d' explosions nucleaires. La representation de la 
difference {mb - Mg) en fonction de nib permet de 
verifier cette hypothese. Neanmoins, le calcul de la 

10 magnitude Ms depend de la periodicite du signal 
enregistre et n'est pas possible en toute rigueur pour 
des evenements regionaux. Au contraire, la presence 
d'une onde de surface, correspondant a une onde de . 
Rayleigh sedimentaire dans des signaux sismiques 

15 proches, caracterise des evenements de nature 
artif icielle. Une methode de detection de ce second 
type d'onde de surface consiste a rechercher 
directement sa presence dans le spectrogramme du 
signal, sa frequence etant connue (entre 0,5 et 1,5 Hz) 

2G et son heure d'arrivee supposee pouvant etre calculee a 
partir de sa vitesse moyenne de propagation et de la 
distance separant 1' epicentre de la station 
d' enregistrement . 

25 Les systemes decrits dans les documents de 

I'art connu ne sont pas operationnellement credibles 
pour plusieurs raisons : 

• Les etudes menees par des geophysiciens , le 
plus souvent riches et detaillees en matiere de 

30 proposition de parametres discriminants, ne proposent 
pas de methode fiable pour 1 ' exploitation automatique 
de ces. parametres . 

• Les etudes menees par des 
informaticiens proposent des systemes qui ne tirent pas- 
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assez parti de la complementarite des donnees et des 
connaissances geophysiques . 

La plupart des etudes de I'art connu 
utilisent des bases de donnees d' evenements sismiques 
5 de taille extremement reduite, ce qui a pour 
consequence de ne pas permettre un apprentissage 
statistique correct. La classification se fait, le plus 
souvent, sur des bases ayant moins d'une centaine 
d ' evenements , comme decrit dans les documents 

10 references [2] et [6] . L'une des bases les plus 
importantes trouvee dans les documents de I'art connu 
est composee de 312 evenements seulement, comme decrit 
dans le document [4] . La consequence directe en est que 
les marges d'erreurs sur les resultats presentes sont 

15 tres elevees, ce qui ne permet pas d'accorder beaucoup 
de confiance a ces resultats. 

La repartition geographique des exemples de 
la base de donnees est un element tres important. La 
plupart des bases regroupent des evenements qui ont eu 

20 lieu dans des regions de taille restreinte (quelques 
dizaines de kilometres de cote), ou les proprietes 
geologiques des sous-sols sont peu diversif iees . La 
recherche de discriminants generaux est done biaisee, 
les discriminants n'etant efficaces que pour une region 

25 donnee. 

De plus, comme decrit dans le document 
reference [1], les evenements des deux classes a 
discriminer peuvent provenir de deux regions 
geographiques clairement distinctes, parfois eloignees 

30 l'une de 1' autre de plusieurs centaines de kilometres. 
II est alors impossible de savoir dans quelle mesure ce 
n'est pas la " coloration " des signaux par les couches 
geologiques traversees qui permet la discrimination, 
plus que les signaux eux-memes . Or les sismologues 

3 5 savent que cette coloration est loin d'etre negligeable 
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et que 1 ' information de localisation est tres 
importante. 

Le plus souvent, un nombre tres limite de 
stations d' enregistrement est utilise. Les signaux sont 
5 . enregistres par deux ou trois stations au maximum, mais 
le plus souvent on se contente d'une unique station. 
L'evenement sismique- est alors represente par un seul 
signal, ce qui reduit considerablement 1 ' information 
utilisable . 

10 Enfin, les evenements integres dans la base 

de donnees sont tres generalement selectionnes selon un 
critere prealablement defini .: magnitude superieure a 
un seuil , rapport signal /bruit superieur a un seuil, 
comme decrit dans le document reference [2]. Mais cette 

15 selection biaise evidemment totalement les resultats. 

Si les etudes de I'art connu ont permis de 
recenser une large gamme de discriminants 
potentiellement utiles a la classification, il s'est 
avere tres difficile de trouver des discriminants 

20 globaux efficaces en raison du nombre eleve de types de 
tirs et de seismes distincts. 

La plupart du temps, les classifieurs 
proposes se revelent inefficaces parce que trop simples 
(separateurs lineaires) ou impossibles a regler du fait 

25 de leur complexite. lis imposent un enorme travail de 
pretraitement sur les donnees, ce qui rend les systemes 
proposes non generalisables . 

Exemple de surveillance sismique 

30 Le laboratoire de detection et de geophysique 

(LDG) du CEA surveille en continu 1 ' activite sismique 
de la terre depuis 1962. Lorsqu'uh evenement sismique 
survient en un endroit du globe, il est enregistre en 
France par un reseau de quarante-deux sismometres 

35 verticaux localises sur le territoire metropolitain, 
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comme illustre sur la figure 1, les stations CP etant 
des stations courtes periodes, et les stations LP des 
stations longues periodes. Une description detaillee du 
reseau de sismometres et de la propagation des ondes 
5 sismiques en France est donnee dans le document 
reference [7]. 

Ce reseau, qui utilisait depuis sa creation 
la transmission par voies hertziennes, est passe 
recemment a une transmission numerique par satellite. 

10 Un filtrage et un gain adaptes du signal permettent de 
detecter les seismes proches ou au contraire les 
seismes lointains de plus longue periode ' appeles 
teleseismes. Le reglage des parametres de filtrage et 
du gain doit permettre de trouver un compromis entre la 

15 detection des evenements sismiques - de magnitude 
relativement faible et le bruit de fond. 

La figure 2 montre des signaux enregistres 
par les sismometres du laboratoire LDG situes entre 84 
et 146 kilometres (soit SBF : 84 km ; PGP : 110 km ; 

20 FRF : 127 km ; LMR : 13 6 km et LRG : 148 km) de 
1' epicentre estime d'un seisme de magnitude 1,9 
localise a 10 kilometres au sud d'Imperia en Italie, le 
9 mai 1996 (heure : 1 heure 0 minute 59 secondes ; 
latitude : 43.34 ; longitude : 8.19 ; magnitude : 1,9). 

25 Sur chaque signal sont repertoriees 

dif f erentes phases sismiques qui serviront lors de 
1 ' analyse detaillee de 1 ' evenement . 

Chaque annee, environ 9000 evenements 
sismiques proches sont ainsi , detectes , dont 800 a 1 200 

30 seismes naturels. 

Les sismologues du laboratoire LDG 
depouillent et analysent quotidiennement les donnees 
enregistrees par les stations du reseau frangais. lis 
publient un bulletin hebdomadaire contenant 1' ensemble 

35 des seismes naturels survenus en France ou dans les 
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regions limitrophes . Un bulletin similaire est publie 
pour les teleseismes. Le tableau I, donne en fin de 
description, est un extrait du bulletin relatif a la 
periode du 9 au 15 septembre 1998. Dans ce tableau on a 
5 les abreviations suivantes : 

HEURE OR : heure d'origine (TU) ; 

LAT : latitude de 1' epicentre (deg) ; 

LON : longitude de 1' epicentre (deg) ; 

PRO : profondeur de 1' epicentre (kxti) ; 

10 ML : magnitude locale ; 

RMS : residu quadratique moyen (s) . 

Ce tableau comporte done 1' ensemble des seismes proches 
detectes durant cette semaine et les caracteristigues 
de chaque evenement : date et heure origine, position 

15 et profondeur de 1' epicentre, magnitude, residu 
quadratique moyen de localisation et region de 
localisation deduite. 

En raison du fort desequilibre entre le 
nombre d'evenements artificiels et celui d'evenements 

20 naturels (en France, les evenements sismiques d'origine 
artificielle sont dix fois plus frequents que les 
evenements sismiques naturels), seuls les evenements 
supposes etre des seismes ou les evenements de elasse 
indeterminee sont extraits du bruit de fond par les 

25 sismologues pour etre ensuite analyses plus precisement 
par un logiciel de localisation. Les autres signaux 
(evenements artificiels pour la plupart) sont archives 
pendant six mois. 

La procedure d' exploitation des signaux 

30 consiste en une localisation deja automatisee, suivie 
d'une phase de caracterisation (detennination de 
1' evenement a I'origine des signaux). 

L' exploitation des ' signaux est effectuee a 
I'aide du systeme global illustre sur la figure 3. 



35 
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L' invention a pour objectif de pallier les 
inconvenients des systemes de I'art connu en proposant 
un nouveau systeme d' intelligence artificielle pour la 
classification d' evenements , d'objets ou de situations 
5 a partir de signaux et de parametres discriminants 
issus de modeles. 

EXPOSE DE INVENTION 

La presente invention concerne un systeme 

10 d' intelligence artificielle pour la classification 
d' evenements, d'objets ou de situations a partir de 
signaux et de parametres discriminants issus de 
modeles, caracterise en ce qu'il comprend au moins une 
branche de traitement comprenant un systeme expert flou 

15 (expert SEF) prenant .une decision a partir de 
proprietes de haut niveau et de parametres 
discriminants de plus bas niveau extraits des signaux 
par des procedures de type traitement de signal, et 
capable d'expliquer sa decision a 1 ' utilisateur par 

20 1 ' intermediaire de regies triees par ordre 
d' applicabilite . 

Dans ce systeme expert flou on effectue une 
descente de gradient sur les parametres : 

• X = y/a 

25 . s = ln(2a^) 

• r = ln(p) 

• d 

avec : 

• y : position des ensembles flous des 
30 premisses 

• a : largeur des ensembles flous des 
premisses 

• p : poids des regies 

• d : degre d' activation de chague classe 
35 pour chaque regie 
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AvcLntageusement le systeme de 1' invention est 
un systeme multi-experts constitue d'au moins deux 
branches de traitement independantes , se configurant 
automat iquement par apprentissage statistique sur des 
5 bases de donnees, ayant des proprietes particulieres, 
et fusionnees par un systeme decisionnel de haut 
niveau, Avantageusement une branche comprend un 
classifieur neuro-flou (expert CNF) prenant ses 
decisions a partir de proprietes de haut niveau et de 

10 parametres discriminants de plus bas niveau extraits 
des signaux par des procedures de type traitement de 
signal. Avantageusement une autre branche comprend un 
reseau neuronal a connexions locales et poids partages 
(expert TDNN) constitues de bancs de filtres adaptatifs 

15 non lineaires, extrayant lui-meme des informations 
discriminantes de representation temps -frequence des 
signaux correspondants a 1 ' evenement . 

L' invention peut etre utilisee dans 
differents domaines d' application et notamment en : 

20 

• SurvBillance evenement s gedphysi<jues 

Le systeme s' applique alors a 1' analyse de 
tout evenement geophysique observable par des signaux 
recueillis par des stations : 
25 - signaux sismiques ; 

- infrasons ; 

- ondes hydro-acoustiques . 

Ces evenements peuvent etre des evenements 
proches (dits regionaux) ou lointains (par exemple- des 
30 teleseismes) , 

la fonction a assurer peut etre : 

- un filtrage pour eliminer les evenements 
non pertinents pour des traitements en aval ; 

- une detection d' evenements particuliers ; 
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- une classification exhaustive en un 
ensemble de groupes d' evenements de meme nature. 

• Surveillance et controle industrials 

5 Le systeme s' applique egalement a 1' analyse 

d'objets ou de processus industriels, pour peu que I'on 
dispose de signaux ou d' images recueillis par des 
capteurs. En voici quelques exemples : 

■ Controle de qualite d'objets ou de produits 
10 manufactures : I'objectif est de controler la forme 

et/ou la position d'objets, detecter et caracteriser 
des defauts. Les experts CNF et SEF utilisent des 
mesures realisees par traitement d' images. L' expert 
TDNN utilise une ou plusieurs images de la piece. 

15 ■ Maintenance predictive d' equipements : 

I'objectif est de prevoir une future def alliance de 
machines , d ' ordinateurs , d ' equipements elec troniques , 
de capteurs afin de donner I'alerte et de pouvoir 
mettre en ceuvre une procedure de correction avant la 

20 panne. Les experts CNF et SEF utilisent des mesures de 
coefficients de haut niveau, des correlations. L' expert 
TDNN utilise des signaux. 

■ Surveillance de processus complexe : 
I'objectif est de surveiller le bon f onctionnement 

25 d'une chalne de production. Les experts CNF et 
SEF utilisent des mesures . de coefficients de haut 
niveau. L' expert TDNN utilise des mesures de 
coefficients de bas niveau. 

30 Dans le domaine geophysique, les proprietes 

de haut niveau peuvent etre la localisation, la 
magnitude, I'heure et la date. Le systeme de 
1' invention permet de realiser la classification 
automatique d' evenements sismiques en trois ensembles : 

35 - seismes naturels ; 
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- explosions (tirs de mines et essais 
militaires) ; 

- coups de terrain (ef f ondrement de couches 

minieres) . 

5 Le systeme s'integre alors dans une chalne de 

traitements automatiques pour operer une fonction de 
filtrage des evenements sismiques. Ses caracteristiques 
principales sont : 

- une fiabilite maximale : le systeme est 
10 capable de prendre des decisions meme avec des donnees 

degradees ou imprecises, voire en 1' absence de 
certaines inf oirmations ; 

- I'acces a 1 ' explication des decisions afin 
de lever un doute eventuel sur une decision. 

15 

BREVE DESCRIPTION DES DESSINS 

La figure 1 illustre le reseau de sismometres 
du laboratoire de detection et de geophysique (LDG) du 
CEA en 1998. 

20 La figure 2 illustre des exemples de signaux 

sismiques enregistres par le reseau de la figure 1. 

La figure 3 illustre un systeme global 
d' exploitation de signaux geophysique de I'art connu. 

La figure 4 illustre le schema de principe du 
25 systeme multi-experts de discrimination d' evenements 
sismiques de 1' invention. 

La figure 5 illustre le schema general 
d'apprentissage du systeme de 1' invention a partir 
d' exemples . 

30 La figure 6 illustre un neurone artificial. 

La figure 7 illustre un reseau de neurones 
artif iciels . 

La figure 8 illustre un classifieur neuro- 

f lou- 
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La figure 9 illustre un mecanisme 
d' activation de cellules de codage. 

La figure 10 illustre des exemples 
d' activation des cellules de codage. 
5 La figure 11 illustre un systeme expert flou 

pour la discrimination d' evenements sismiques. 

Les figures 12A a 12D illustrent des 
pretraitements successifs appliques a un signal 
sismique . 

10 La figure 13 est une architecture d'un reseau 

de neurones a connexions locales et poids partages. 

La figure 14 illustre la representation des 
epicentres des evenements sismiques de 1962 a 1995. 

EXPOSE DETAILLE DE MODES DE REALISATION 
Description generale du systeme 

Le systeme de I 'invention comporte au moins 
une branche de traitement contenant un systeme expert 
flou. Lorsqu'il comprend plusieurs branches de 
traitements independantes , on parle de systeme multi- 
experts . 



15 



20 



• Prise de decision mult i -experts 

Le principe de la prise de decision multi- 
25 experts, qui est done I'un des modes de realisation de 
1' invention, est 1 ' exploitation de la synergie entre 
plusieurs branches de traitement complementaires . Cette 
complementarite reside dans : 

A. Les performances generales : une branche 
30 est plutot generaliste (performances assez 

bonnes dans la plupart des cas), une autre 
est plutdt specialiste (tres bonnes 
performances sur certains cas difficiles, 
taux d'erreurs plus eleve sur les cas ne 
35 relevant pas de sa competence) . 
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B. Les performances selon le cas traite : urie 
branche peut gtre plus apte qu'une autre a 
traiter tel cas particulier. 

C. La nature des entrees (signaux ou donnees 
5 de haut niveau) . 

D. La nature des sorties (simple donnee de la 
classe, estimation de la certitude de la 
decision, explication formelle de la 
decision) . 

10 La figure 4 illustre le schema de principe du 

systeme multi-experts de discrimination d'evenements 
sismiques de 1' invention. Ce systeme est constitue de 
plusieurs branches de traitement independantes , chacune 
ayant des proprietes particulieres , fusionnees par un 

15 systeme decisionnel de haut niveau. 

Ces branches sont : 

~ un classifieur neuro-flou, dit CNF, prenant 
ses decisions a partir de proprietes de haut niveau des 
evenements (par exemple pour les evenements sismiques : 
20 localisation, magnitude, heure, jour de la semaine) et 
des parametres de plus has niveau extraits des signaux 
par des procedures de type traitement du signal ; 

- un systeme expert flou, dit SEF, prenant 
une decision de fagon independante a partir des memes 

25 informations, et capable d'expliquer sa decision a 
1 'utilisateur par 1 ' intermediaire de regies triees par 
ordre d' applicabilite a I'evenement traite ; 

- un reseau neuronal a connexions locales et 
poids partages, dit TDNN, constitue de bancs de f litres 

30 adaptatifs non lineaires, extrayant lui-meme des 
informations pertinentes de representations temps- 
frequence des signaux correspondant a I'evenement. 

Ces trois branches se configurent 
automatiquement par apprentissage statistique sur des 

35 bases de donnees d'evenements sismiques. 
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• Apprentissage a partir examples 

L'apprentissage a partir d'exemples consiste 
a construire un modele du systeme de prise de decision 
5 par ajustement progressif de parametres a partir de 
donnees. Ce modele doit etre capable d'associer la 
bonne decision (sortie) a un ensemble de donnees 
decrivant le cas traite (entrees) . Cela se fait 
progressivement , par presentation iterative des cas 
10 disponibles dans la base d'exemples en" entree du 
systeme. Une telle procedure est illustre sur 
1 ' organigramme de la figure 5. 

Dans 1' invention, le modele qui apprend peut 
etre soit un reseau de neurones artificiels, soit un 
15 systeme expert flou. 

Une fois que le systeme a termine sa phase 
d' apprentissage, ses parametres internes sont figes et 
le systeme est pret a etre utilise. 

20 • Reseau de neurones de type « Perceptron multicouche » 

Un tel reseau de neurones artificiels de type 
« Perceptron multicouche » est un modele particulier de 
reseau de neurones pouvant etre utilise comme systeme 
de prise de decision. II est constitue d'un reseau 
25 d' automates de calcul simples, les " neurones 
artificiels 

Un neurone N"', comme illustre sur la figure 

6, est une entite constitute d'un vecteur-poids 
Wy={vv,y} et d'une fonction de transfert non lineaire (|). 

30 II admet en entree un vecteur X = {x,} et effectue une 

f \ 



transformation de ces entrees de type yj^(p 



Par similitude avec, le vocabulaire employe en 
neurophysiologie, on dit que chaque entree Xi est 
reliee au neurone N"^ par une connexion synaptique. Un 
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poids synaptique Wij module I'efficacit^ de cette 
connexion. 

Dans un reseau de neurones artificiels, comme 
illustre sur la figure 7 , les neurones sont assembles 
5 en couches successives. Une couche est definie comme un 
ensemble de neurones n'ayant pas de connexions entre 
eux, mais pouvant avoir des connexions avec des 
neurones des couches precedentes (entrees) ou suivantes 
(sorties) . En general, on ne connecte que des neurones 
10 de couches successives. 

L' apprentissage consists a modifier 
progressivement les valeurs des poids Wij jusqu'a ce que 
les sorties du reseau, qui est constitue d'un certain 
nombre de couches de neurones, correspondent aux 
15 sorties desirees. 

Pour cela, on definit une erreur de 
classification que 1 ' on cherche a minimiser. L' erreur 
la plus couramment utilisee est 1 ' erreur quadratique, 

definie par £= V v-^/t'^^A- ^^^^^^ f . La methode consists a 

2 0 faire une descents de gradient sur les poids par 

dE 

I'equation Aw- =-a avec a> 0. Cette equation, 

lorsqu'on la developpe, donne la formule de correction 
pour chaque poids du reseau. 

25 • Systeme expert flou 

Un systeme expert flou est un autre models ds 
systeme de prise de decision. II pressnte I'avantags 
sur Is rsseau ds neurones de donner une forme 
d' explication de ses decisions. II est constitue d'un 
30 ensemble d'unites de calcul, les regies d' inference 
f loue " . 

Une regie d' inference floue est une entite ds 
la forme si < premisse > alors < conclusion > La 
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premisse est la partie que I'on cherche k mettre en 
correspondance avec les donnees d' entrees. 

En effet, un ensemble flou est un ensemble 
dont les frontieres sont progressives, au contraire 
5 d'un ensemble classique, dont les frontieres sont 
nettes. Un element fait done plus ou moins partie de 
chaque ensemble flou. Lorsque les donnees sont de 
dimension 1 , un ensemble classique peut etre represente 
par un rectangle (appartenance = 1 a I'interieur, 0 au- 

10 dehors), alors qu'un ensemble flou peut etre un 
triangle, un trapeze, une gaussienne... 

De meme que precedemment , 1 ' apprentissage 
consiste a modifier progressivement les valeurs des 
parametres jusqu'a ce que les sorties du systeme expert 

15 flou correspondent aux sorties desirees. 

Quatre types de parametres sont calcules par 
apprentissage : la position et-la largeur des ensembles 
flous des premisses, les poids des regies et le degre 
d' activation de chaque classe pour chaque regie. 

20 En phase d' utilisation operationnelle , le 

systeme expert flou fournit, en plus de la classe 
attribute a 1 ' evenement sismique , la liste des regies 
applicables par ordre decroissant de pertinence. 
Certaines de ces regies peuvent etre en contradiction 

25 avec les autres, ce qui permet d' examiner les solutions 
alternatives, mais c'est I'agregation du resultat de 
toutes les regies qui fournit le resultat global. 

L' utilisateur a done a sa disposition : 

- la decision brute (seisme, explosion, coup 
30 de terrain) ; 

- la liste des regies applicables ; 

- la liste des regies apportant une 
contradiction a cette decision ; 
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- la raison de la decision de chaque regie 
(par examen de la coherence entre les donnees et les 
ensembles flous correspondants) . 

Un exemple de regie de decision trouvee par 
5 le systeme est le suivant : 

si (Heure est milieu jde apres-midi) 

et (Latitude est tres ^)roche_de 43.5° N) 

et (Longitude est tres _j)roche_de 5.5° E) 

et (Magnitude est environ 2.7) 

10 et (Date est ^/e_/?re/i?>^/2ce samedi) 

alors (avec niveau de confiance = 0,8) 

(seisme est improbable) 
(explosion est probable) 
(coup de terrain est improbable) 

15 

Dans 1' invention, pour des raisons de 
convergence difficile, on parametre cette descente de 
gradient en introduisant des variables intermediaires . 
Si I'on veut effectuer une descente de gradient sur un 

20 parametre p avec p=<t){s), ' ^ etant une fonction 
derivable, strictement monotone, independante de p et 
des valeurs des exemples servant a 1 ' apprentissage , on 
a les memes solutions finales en effectuant une 
descente de gradient sur s. L'avantage d'un tel 

25 changement de variable est qu'il devient possible de 
changer la fagon d'atteindre la solution, et en 
particulier de faciliter la convergence dans les cas 
dif f iciles . 

Dans 1' invention, on optimise les parametres • 

30 suivants : 

(1) La position y des ensembles flous des 
premisses : lorsque la descente de gradient est 
appliquee directement sur ce parametre, on obtient 
generalement une convergence difficile. Ceci s'explique 
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par le fait que la variation de la position y des 
ensembles flous des premisses n'est pas une fonction 
croissante de la distance a I'exemple. On corrige ce 
phenomene en posant x=y/a. 
5 (2) La largeur a des ensembles flous des 

premisses : lorsque les donnees sont structurees en 
groupes de tallies tres differentes, I'algorithme ne 
peut pas converger. En etudiant la variation relative 
Aa/a, on decouvre qu'elle n'est pas bornee (c'est-a- 

10 dire que rien ne I'empeche de tendre vers des valeurs 
infinies) . Lorsque les donnees sont tres groupees, 
cette variation prend ef f ectivement des valeurs tres 
elevees . Pour avoir une modification relative plus 
faible lorsque les donnees sont rapprochees, on pose 

15 s=ln(2a^). 

(3) Les poids p des regies : c'est le 
parametre le plus difficile a regler. Avec une descente 
de gradient directs, les poids les plus faibles 
diminuent et deviennent negatifs, ce qui leur fait 

20 perdre toute signification et fait diverger 
I'algorithme. On choisit done une fonction d' activation 
positive en imposant une contrainte supplementaire : 
pour differents examples avec le meme niveau activation 
de la regie, la variation de ce niveau doit etre la 

25 meme si les conclusions sont egales . La consequence est 
que la variation relative des poids des regies doit 
etre constante lorsque les exemples ont le meme degre 
d' appartenance aux ensembles flous. Ceci est realise en 
posant r=ln{p) . 

3 0 (4) Le degre d' activation d de chaque classe 

pour chaque regie. 

La descente de gradient est done effectuee 
non pas sur y, a et p, mais sur : 
■ x=y/a 
35 ■ s=ln{2o^) 
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• r=ln(p) 

Pour d, on n'effectue pas de changement de variable. 
Ces changements de variables assurent une tres bonne 
gualite de convergence et permettent I'obtention de 
5 systemes -experts flous tres efficaces. 

• Base examples et validation 

La base d'exemples utilisee doit verifier 
deux principes f ondamentaux : 

10 - etre gualitativement representative du 

probleme reel (distribution des exemples conformes a la 
distribution reelle) ; 

- etre quant itacivement representative du 
probleme (nombre d' exemples suffisant pour constituer 

15 un echantillonnage satisf aisant ) . 

II existe plusieurs methodologies 

d' apprentissage et de validation. Dans la procedure la 
plus simple, on divise la base d' exemples en deux bases 
disjointes : la base d' apprentissage et la base de 

20 test. . On entraine le systeme par apprentissage sur la 
premiere et on verifie son bon , f onctipnnement sur la 
seconde. Une base d' exemples ne verifiant pas 
conj ointement les deux proprietes enoncees ci-dessus 
fait courir le risque de conduire a un systeme 

25 incapable de generaliser correctement , c'est-a-dire de 
fonctionner sur des exemples nouveaux, non presentes au 
cours de la phase d' apprentissage , 

Que ce soit pour les reseaux de neurones 
artificiels, pour les systemes-experts flous ou plus 

30 generalement pour tout systeme congu par apprentissage 
statistique sur des donnees experimentales , il est 
primordial d'utiliser une base d' exemples suffisante en 
qualite et en quantite, et de valider le systeme mis au 
point par des procedures rigoureuses . 



35 



wo 01/03064 



PCT/FROO/01804 



21 

Description detaillee 

• Description d'un evenement sismi<^e 

Un evenement sismique a identifier peut etre 
decrit par : 

5 . - la donnee des signaux en provenance du 

reseau de stations sismiques, ou 

- des proprietes de plus haut niveau, 
directement mesurables ou calculees par des modeles 
geophysiques . On peut utiliser ^ par exemple la 

10 localisation de 1 ' evenement (latitude et longitude), sa 
magnitude et le moment ou il est survenu (heure et jour 
de la semaine) . Par exemple les informations de haut 
niveau sont les suivantes : le jeudi 7 avril 1966 a 12 
heures, un seisme de magnitude 1,4 s'est produit par 

15 02^35' 06" de longitude est et 49<=*12'25" de latitude 
nord. 

• Le classifieur neuro-flou 

Le classifieur neuro-flou (expert CNF) , comme 

20 illustre sur la figure 8, est constitue d'un codage 
neuro-flou des donnees suivi d'un perceptron 
multicouche. II s' applique aux donnees de haut niveau. 

Le codage neuro-flou consiste a associer 
plusieurs cellules de codage a chaque variable (ou 

25 groupe de variables) d' entree, chaque cellule ayant une 
region d' influence modelisee par une fonction 
definissant son mecanisme d' activation . La presentation 
d'un vecteur de valeurs genere alors un schema 
deactivation des cellules de codage associees. 

3 0 La figure 9 illustre ce mecanisme 

d' activation des cellules de codage. La presentation 
d'une valeur genere un schema d' activation 
correspondant a la reponse ^ impulsionnelle de chaque 
fonction d' activation a la valeur presentee. Les 

35 niveaux de gris affectes aux centres des cellules 
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indiquent leur niveau deactivation compris entre 0 et 1 * 
(noir : 1, blanc : 0) . 

• La figure 10 montre des exemples de schemas 
deactivation generes par la presentation de valeurs 
5 typiques. II s'agit d'un codage de type curseur. La 
valeurs faibles (respectivement elevees) activent 
pref erentiellement les cellules de gauche 
(respectivement droite) . 

L'interet de ce codage est multiple : 
10 - par sa nature meme, il permet de 

representer des donnees incompletes, imprecises ou 
incertaines at de les utiliser efficacement pour la 
prise de decision ; 

- par ses proprietes de conditionnement non 
15 lineaire des donnees, il facilite les traitements 
ulterieurs (ici, la classification} . 

Ce codage neuro-flou est realise en plusieurs 
etapes successives : definition des sous-groupes de 
caracteristiques , choix et placement des cellules de 
20 codage affectees a chaque groupe, determination des 
parametres de la region d' influence de chaque cellule. 
Les details de cette procedure sont exposes dans le 
document reference [8] . 

Une fois les donnees codees, elles sont 
25 analysees par le perceptron multicouche qui en calcule 
la classe. 

• Le systems expert flou 

Dans un mode de realisation le systeme de 
30 . 1' invention comporte une seule branche de traitement 
basee sur un tel systeme expert flou 

Le systeme expert flou (expert SEF) 
s' applique egalement aux donnees de haut niveau. 

Dans la figure 11, on represente un systeme 
35 expert flou a cinq regies (une regie par ligne) . 
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Pour chaque ligne, les cinq, colonnes de 
gauche representent les premisses et" correspondent a 
cinq entrees : I'heure, la latitude, la longitude, la 
magnitude et la date. Les premisses sont composees 
5 d' ensembles flous en forme de gaussiennes qui 
recouvrent les domaines de la variable d' entree 
conduisant a un renforcement de I'activite de la regie. 

Quatre types de parametres sont calcules par 
apprentissage : la position et la largeur des ensembles 
10 flous des premisses (colonnes 1 a 5) , les poids des 
regies (colonne 6), qui permet de specifier le degre 
d' importance de chaque regie dans le processus de 
decision et le degre deactivation de chaque classe 
(seisme naturel, explosion ou coup de terrain) pour 
15 chaque regie (colonnes 7 a 9) . 

A chaque presentation d'un exemple a 
classifier, on calcule le contenu de la colonne 10 et 
de la colonne 6 : 

■ La colonne 10 donne 1' activation de chaque 
20 regie (et permet done d'estimer son adequation avec le 

cas traite) . 

■ La ligne 6. est la synthese des decisions 
des cinq regies et donne la reponse globale du systeme 
expert flou (ici, la decision est « explosion ») . cette 

25 synthese est faite en calculant le barycentre des 
decisions de toutes les regies (colonnes 7 a 9) 
ponderees par le niveau d' activation correspondant 
(colonne 10). Sur la figure 11, .la position du 
barycentre pour . chaque classe est symbolisee par un 
30 trait vertical ligne 6, colonnes 7^9. 

L' apprentissage est realise en deux temps : 
- une premiere phase consiste a positionner 
les ensembles flous (centres et largeurs), par exemple 

r 

au moyen d'un algorithme dit des C-moyennes floues, 
35 comme decrit dans . le document reference [9] ; 
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- une deuxieme phase consiste a realiser une 
descente de gradient sur les quatre types de 
parametres. 

5 • Ire reseau neuronal a connexions locales et poids 
partages 

Au contraire des deux branches precedentes, 
le reseau de neurones a connexion locale et poids 
partages (expert TDNN) admet en' entree les signaux 

10 sismiques eux-memes et apprend a extraire par 
apprentissage non seulement la procedure de decision, 
mais aussi les parametres discriminants qui vont servir 
de base a cette decision. Ce reseau neuronal est de 
type perceptron multicouche a connexions locales et 

15 poids partages prenant en entree les spectrogrammes 
pretraites des signaux sismiques, comme decrit dans le 
document reference [10], Ces . spectrogrammes sont 
obtenus. par 1 ' application d'une transformee de Fourier 
a fenetre glissante sur le signal. 

20 Les figures 12A a 12D illustrent les 

pretraitements successifs appliques a chaque signal 
sismique aboutissant a un spectrogramme final a 15 
bandes f requentielles : la figure 12A illustre le 
signal initial ; la figure 12B illustre le 

25 spectrogramme deduit du signal avec 50 bandes 
f requentielles ; la figure 12C illustre le 
spectrogramme « debruite » ; la figure 12D illustre le 
spectrogramme apres reduction de 50 a 15 bandes 
f requentielles . 

30 Le spectrogramme obtenu. est ensuite pretraite 

puis presente en entree d'un reseau de neurones de type 
TDNN. Chaque reseau est specialise dans le traitement 
des signaux enregistres par une station donnee. 

La figure 13 illustre 1 ' architecture d'un 

35 reseau TDNN specialise dans la classification des 
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spectrograiranes deduits des signaux enregistres par un 
sismometre donne, ce reseau comportant quatre couches 
de neurones'. La couche d' entree a des connexions 
locales et des poids partages (4 trames avec un 
5 decalage de 2 trames) avec la premiere couche cachee. 
Celle-ci a de meme des connexions locales a poids 
partages (9 trames avec un decalage de 5 trames) avec 
la seconde couche cachee, totalement connectee a la 
derniere couche. 

10 Les poids partages rendent 1 ' architecture 

plus robuste aux petits decalages des pointages des 
phases ou de trames manquantes ou erronees . Cependant, 
a cause de la vitesse de propagation des ondes P 
(compression) et S (cisaillement ) , la duree entre 

15 I'arrivee de la phase P et de la phase S varie en 
fonction de la distance entre la station 
d' enregistrement et 1' epicentre de I'evenement, ce qui 
complique 1 ' apprentissage . La solution adoptee consiste 
a aligner le pointage de la phase P sur la 10^""® trame 

20 et celui de la phase S sur la 60^""^ trame. 

• La prise de decision finale 

Pour la prise de decision finale, on suppose 
que toutes les sorties sont comprises dans I'intervalle 

25 reel[-M]. Cette prise de decision consiste en une 
association des reponses fournies par les trois 
branches afin d'augmenter la fiabilite, Elle peut etre 
realisee pair un calcul de moyennes arithm^tigues sur 
les sorties homologues de chacune des trois branches. 

30 Pour chacune des trois sorties Si du systeme global, on 
a alors : 

i=1...3 

La certitude sur la reponse est evaluee par 
un coefficient calculable seulement si le systeme est 
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en situation de prise de decision (c ' est-a-dire s'il y 
a une et une seule sortie strictement positive) . Ce 
coefficient est alors egal a^ la moyenne des valeurs 
absolues des sorties : 




K<0.2 : doute 

/Te] 0.2, 0.4] : prudence 

ATejOAO.d] : certitude raisonnable 

^6]o.6,0.8] : grande certitude 

10 • K>O.E : certitude quasi absolue 



Ainsi, on peut avoir par exemple : 



Systeme 


Ciasse choisie 


Detail des reponses par classe 


Degre de ceniiude 


1 


Classe 3 


(-0,9 -0.4 +0,8) 


0,7 : grande certitude 


2 


Ind^termination 1 ou 3 


(+0,1 -0,5 +0,3) 


Incertitude complete 


3 


Classe 1 


(+0,2-0,6 -0.2) 


0,3 : prudence 


Fusion 


Classe 3 


(-0.2 -0,5 +0,3) 


0,3 : prudence 



15 

Exemple de mise en oeuvre de 1 ' invention : 
discrimination d ^ evenements sismiques regionaux 
• Localisation de I 'evenement 

La discrimination « evenement naturel / 

20 evenement artificiel » . est une etape majeure de la 
surveillance sismique, realisee rapidement des la 
lecture des signaux lors de la phase de depouillement, 
puis affinee au cours de chaque nouveau traitement. On 
evalue a environ sept ans le temps necessaire aux 

25 analystes pour §tre veritablement operationnels . Ces 
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analystes etant de vsritables experts, il est difficile 
d'exprimer clairement leur methode de raisonnement 
basee a la fois sur des connaissances expertes et sur 
un raisonnement par cas. 
5 La localisation d'un evenement sisinique est 

obtenue apres une succession d' aller-retours entre deux 
phases principales : 

- le pointage des differentes phases 
sismiques effectue sur les signaux enregistres par les 

10 stations ayant detecte 1' evenement et le calcul de la 
magnitude ; 

- la localisation elle-m§me, realisee a 
I'aide d'un modele mathematique congu par les 
sismologues. 

15 A 1' issue de la phase de pointage, la 

localisation de 1' evenement peut-etre effectuee a 
I'aide d'un logiciel de simulation connu de I'homme de 
metier. II fait appel a des modeles sismologiques 
contenant des informations sur la vitesse des ondes, 

20 les differents types d'ondes et leur mode de 
propagation, des corrections d' altitude selon les 
stations... Plusieurs hypotheseis de localisation sont 
proposees, associees a un degre de coherence des 
donnees. Si celles-ci ne conviennent pas a 1' expert, il 

25 modifie ses pointes de phases puis relance une 
localisation. Ce cycle est reitere jusqu'a I'obtention 
d'un resultat juge satisf aisant . 

La qualite de la localisation depend du 
nombre et de la qualite des stations ayant servi a la 

30 localisation de 1' evenement, ainsi que de leur 
repartition azimutale. Les evenements situes en France 
sont done, en regie generale, mieux localises que les 
evenements etrangers . Pour af finer la localisation dans 
ce dernier cas,^ les experts accedent regulierement a 

35 des donnees etrangeres. Avec les seules donnees 
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frangaises, la precision de localisation des evenements 
survenant en France est en moyenne de cinq kilometres. 
Elle est estimee a environ un kilometre dans le 
meilleur des cas. 
5 Le tableau II en fin de description donne un 

exemple de resultats fournis par la procedure de 
localisation. La.partie superieure resume les resultats 
obtenus : a chaque valeur d'heure origine, de magnitude 
et de localisation (latitude-longitude) est associee le 

10 niveau d' imprecision estime. La partie inferieure 
indique les stations ayant servi a la localisation et 
les erreurs moyennes quadratiques (RMS) obtenues en 
fonction de 1' hypo these de profondeur de 1' epicentre. 
Dans le cas des coups de terrain, la profondeur est 

15 arbitrairement fixee a un kilometre. 



• Caracterisation de 1 ^ evenement 

On s'interesse ici aux evenements sismiques 
locaux et regionaux, c'est-a-dire survenus en France 
20 metropolitaine ou dans les regions limitrophes . Ces 
. evenements sont souvent qualifies de proches par 
opposition aux teleseismes lointains (epicentre situe a 
plusieurs milliers de kilometres du capteur) . 

Trois types d' evenements sismiques doivent 
25 etre discrimines : 

- les seismes, evenements sismiques d'origine 
naturelle ; 

- les explosions terrestres (tirs de mines, 
de carrieres, de chantiers...) ou en mer (deminage, 

3 0 essais d'armes...) ; 

- les coups de terrains correspondant a 
I'effondrement d'une couche miniere et lies a 
1' exploitation de la mine. 

l: analyse de I'etat de I'art a montre I'echec 
35 des approches basees sur une discrimination a partir de 
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signaux sismiques uniquement. On exploite done 
1' ensemble des donnees disponibles en adoptant une 
approche basee sur la fusion multi -experts et multi- 
sources. La conception du systeme de discrimination 
• 5 automatigue d'evenements sismiques est basee sur trois 
modules : 

- les deux premiers (experts CNF et SEF) 
sont des modules effectuant une discrimination a partir 
des seules donnees de haut niveau, deduites par le 

10 modele inverse du laboratoire LDG. Aucun signal 
sismique n'est done directement pris en compte a ce 
niveau ; 

- le troisieme (expert TDNN) est base sur 
1' analyse des signaux sismiques. 

15 

• Donnees utillsees 

Repartition geographique des evenements 

Les evenements sismiques a analyser se 

repartissent sur 1' ensemble du territoire frangais et 
20 dans un perimetre proche. Les epicentres des evenements 

enregistres par le laboratoire LDG. de 1962 a 1996 sont 

illustres sur la figure 14. 

Donnees de haut niveau 

25 Chaque evenement sismique est caracterise 

par les informations suivantes : la date et I'heure 
origine de 1' evenement, la latitude et la longitude de 
1' epicentre, et sa magnitude. 

L'heure et la date sont retenues en raison 

30 des regies regissant les tirs de carrieres ou de mines 
en France, interdisant les tirs nocturnes ou durant les 
week-ends et jours feries. Neanmoins, des derogations 
sont accordees a certains chantiers, par exemple pour 
eviter les nuisances a la circulation. 
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La magnitude est retenue car les coups de 
terrain ont une magnitude typique d'apres les 
sismologues (environ 3). De plus, seuls les seismes 
peuvent fournir les plus fortes magnitudes. Plusieurs 
5 valeurs de magnitude sont prises en compte lorsqu'elles 
sont disponibles . 

La localisation de 1' epicentre, caracterisee 
par sa latitude et sa longitude, est egalement une 
caracteristique majeure. II existe neanmoins certaines 
10 mines situees dans des regions a forte sismicite et 
susceptibles d' occasionner des coups de terrain. 

Donnees de bas niveau 
Les donnees de bas niveau sont les signaux en 

15 provenance des 42 stations sismiques du laboratoire LDG 
(voir la figure 1) . Les pretraitements portent 
essentiellement sur la creation des spectrogrammes des 
signaux sismiques, qui sont non stationnaires . Ces 
spectrogrammes sont obtenus par 1 ' application d'une 

20 transformee de Fourier a fenetre glissante sur le 
signal. Dans un premier temps, le signal, echantillonne 
a 50 Hz , ■ est segmente en trames de deux secondes 
decalees d'une seconde par une fenetre de Hamming. On 
calcule ensuite la densite spectrale d'energie sur 50 

25 bandes f requentielles avec elimination de la composante 
continue. Puis, on applique une transformation 
logarithmique avec soustraction de bruit, dont le 
modele est suppose logarithmique, dans chaque 
bande selon 1' equation max(ln(l + ;c)-/i(&nm)-o'(/7rmr),0) , ou 

30 ^ifbruit; et a(bruit) correspondent a la moyenne et a 
I'ecart-type du bruit estimes sur une periode 
anterieure au pointage de I'onde P. Enfin, on reduit le 
nombre de bandes f requentielles de 50 a 15 par une 
compression pseudo-logarithmique des hautes frequences. 



35 
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• Resultats obtenus 

Le systeme decrit ci-dessus classifie les 
5 evenements sisitiiques regionaux frangais avec les 
performances suivantes : 

- 86 % pour les seismes ; 

- 91 % pour les explosions ; 

- plus de 99% pour le coups de terrain. 
10 Les performances globales sont d' environ 90%. 
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REVENDICATIONS 

1. Systeme d' intelligence artificielle de 
classification d' evenements, d'objets ou de situations 

5 a partir de signaux et de param^tres discriminants 
issus de modeles, caracterise en ce qu'il comprend au 
moins une branche de traitement comprenant un systeme 
expert flou prenant une decision a partir de proprietes 
de haut niveau et de parametres discriminants de plus 
10 has niveau extraits des signaux par des procedures de 
type traitement de signal, et capable d'expliquer sa 
decision a 1 'utilisateur par 1 ' intermediaire de regies 
triees par ordre d' applicabilite . 

2. Systeme salon la revendication 1, dans 
15 lequel, dans le systeme expert flou on effectue une 

descents de gradient sur les parametres : 

• X = y/a 

• s = ln(2G^) 

• r = ln{p) 
20 • d 

avec : 

• y : position des ensembles flous des 
premisses 

a : largeur des ensembles flous des 
25 premisses 

p : poids des regies 

d : degre d' activation de chague classe 
pour chaque regie 

3. Systeme selon la revendication 1, qui est 
3 0 un systeme multi-experts constitue d'au moins deux 

branches de traitement independantes, se configurant 
automatiquement par apprentissage statistique sur des 
bases de donnees, ayant des proprietes particulieres et 
fusionnees par un systeme decisionnel de haut niveau. 
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4. Systeme selon la revendication 3, dans 
lequel une branche comprend un classifieur neuro-flou 
prenant ses decisions a partir de proprietes de haut 
niveau et de parametres discriminants de plus bas 

5 niveau extraits des signaux par des procedures de type 
trait ement de signal. 

5. Systeme selon I'une quelconque des 
revendications 3 ou 4, dans lequel une branche comprend 
un reseau neuronal a connexions locales et poids 

10 partages constitues de bancs de filtres adaptatifs non 
lineaires, extrayant lui-meme des informations 
. discriminantes de representation temps -frequence des 
signaux correspondants . 

6 . Systeme selon 1 ' une quelconque des 
15 revendications precedentes, qui est un systeme de 

classification d'evenements geophysiques . 

7. Systeme selon la revendication 6, dans 
lequel les proprietes de haut niveau sont la 
localisation, la magnitude, I'heure et la date. 



wo 01/03064 



PCT/FROO/01804 




wo 01/03064 



PCT/FROO/01804 



Osl 

eg 

Li. 




-^/12 




wo 01/03064 



PCT/FROO/01804 




wo 01/03064 



PCT/FROO/01804 



par 


erts 


o 


Q. 


bal 


-ex 


o 










E 


o 






c 


'cn 


o 




'So 





















CD 


CO 


c 






es 






o 


i2 








1 
















CD 


<o 

CD 


/Ve 




c: 




««»^ 


c 


:3 


o 


CD 







CD 



wo 01/03064 



PCT/FROO/01804 



c 



Definition de la structure du systeme 



^Initiali 



Initialisation des parametres internes du systeme 



c 



Presentation d'un exemple de la base 



^ Calcul de la sortie (decision^ 



( 



Correction des parametres 
internes du systeme 
% 



Calcul de Terreur par rapport a la sortie desiree 



esiree^ 




Fin de Tapprentissage 



non 



FIG. 5 




FIG. 6 







— ► 







wo 01/03064 PCT/FROO/01804 




wo 01/03064 



PCT/FROO/01804 



7/1-2 




wo 01/03064 



PCT/FROO/01804 



60 

O 
O 

0) 
0) 



o 



> 
cs 



1 

0,9 
0,8 
0,7 
0,6 
0,5 
0,4 
0,3 
0,2 

0,1 
0 



-©- 



'2 X "0 2 ' . 4 
valeur de la caracteristique 



FIG. 9 



valeur faible 
valeur moyenne 
valeur elevee 



© 


0 


0 


0 


0 0 


0 


0 


• 


0 


0 0 


0 


0 


0 


0 


• 0 



FIG. 10 



wo 01/03064 



PCT/FROO/01804 




wo 01/03064 



11/1:2 



PCT/FROO/01804 



in 

s 



03 

c/3 D 



rr) <D 



1 I 



CO 
CD 




0) 

=3 



00 



00 



sduisj 



wo 01/03064 



PCT/FROO/01804 



1^ 



CD 




apnjijDj 



INTERNA! lAL SEARCH REPORT 



Inten tu,, Applleatlon No 

PCT/FR 00/01804 



A. CLASSinCATlON OF SUBJECT MATTER 

IPC 7 G06N5/04 



Accorxjrng to Intemationai Patent Classification (IPC) or to both national dassification and IPC 



B. FIELDS SEARCHED 



Minimum documentation searched (classification system followed by classification symbols) 

IPC 7 606F G06N GOIV 



Documentation searched other than minimum documentation to the extent that such documents are included in the fieids searched 



Electronic data base consulted during the international search <name of data base and, where practical, search terms used) 

EPO-Internal 



C. DOCUMENTS CONSIDERED TO BE RELEVANT 



Category * Citation of document, with indication, where appropriate, of the relevant passages 



Relevant to claim No. 



CHU CH -K P ET AL: "FIRST BREAK 
REFRACTION EVENT PICKING USING FUZZY LOGIC 
SYSTEMS" 

IEEE TRANSACTIONS ON FUZZY SYSTEMS, US, IEEE 
INC. NEW YORK, 
vol . 2, no. 4, 

1 November 1994 (1994-11-01), pages 

255-266, XP000469682 
ISSN: 1063-6706 
abstract 

page 257, left-hand column, line 45 -page 

258, right-hand column, line 36 

page 265, left-hand column, line 15 - line 

37 

EP 0 468 229 A (HNC INC) 

29 January 1992 (1992-01-29) 

column 13, line 15 -column 15, line 9 

-/-- 



1,4.6 



1,4 



Further documents are listed in the continuation of box C. 



Patent fajnily members are listed in annex. 



* Special categories of cited documents : 

'A' documerrt defining the general state of the art which is not 

considered to be of paiticular relevance 
*E' earlier document but putiTished on or after the intemationai 

filing date 

"L" document wtTich may throw doutrts on priority ciaim(s)or 
which is cited to establish the publication date of another 
citation or other special reason (as specified) 

*0' documertt referring to an oral disclosure, use. exhibition or 
other means 

"P' document pubfished prior to the international fSing date but 
later ttwi the priority date claimed 



f later document pulDlished after the intemationai filing date 
or priority date and not in conflict with the application but 
cited to understand ttie principle or ttieoiy urtdeitying the 
invention 

"X" document of particular relevance: the claimed invention 
cannot be considered novel or cannot be considered to 
involve an Inventive step wfien the document is tal<en alone 

*Y' document of particular relevance; the claimed invention 

cannot be considered to involve an inventive step when the 
document is combined with one or more other such docu- 
ments, such combination being obvious to a person skilled 
in the art 

*&' document member of the same patent family 



Date of ttie actual completion of ihe intematiorial search 



5 September 2000 



Date of mailing of the interriational search report 



12/09/2000 



Name and mailing address of the ISA 

European Patent Office, P.B. 5818 Patentlaan 2 
NL - 2280 HV Rijswijk 
Tel. (431-70) 340-2040. Tx. 31 651 epo nl. 
Fax: (+01-70)340-3016 



Authorized officer 



Kingma, Y 



Form PCT/1SA«10 (saoorvl sheet) (July 1992) 



INTl ATIONAL SEARCH REPORT 



2 



1 IrHBTT -lal MpfNIGBllUil IW 

i DTT /ro nA/nioA>i 
1 rUl/rK UU/UIoU4 


C.(Contlnuatlon) DOCUMENTS CONSIDERED TO BE RELEVANT 


Categoiy • 


Citation of document, with indication.wtiere appropriate, pi the relevant passages 


Relevant to daim No. 


A 


MAGOTRA N ET AL: "SEISMIC SIGNAL 

DISCRIMINATION USING ADAPTIVE SYSTEM 

PARAMETERS" , PROCEEDINGS OF THE MIDWEST 

SYMPOSIUM ON CIRCUITS AND SYSTEMS, US, NEW 

YORK, IEEE, VOL. CONF. 33, PAGE(S) 84-87 

XP000280121 

the whole document 


1,4-6 



Form PCT/ISA«10 (oontinualion of second stieet) Wulv 1 992> 



INTERNATH iL SEARCH REPORT 

information on patent family membefs 



Intent nai Appllcatton No 

PCT/FR 00/01804 



Patent document 
cited in search report 



Publication 
date 



Patent family 
me(nber(s) 



Piiblicalion 
date 



EP 0468229 A 



29-01-1992 



JP 
US 



5151188 A 
5398300 A 



18-06-1993 
14-03-1995 



Fonn PCT/ISA«10 (patent tamity ainax) (Jiiy 1992) 



RAPPORT DE RECHE HE INTERNATIONALE 



Denu Intematlonale No 

PCT/FR 00/01804 



A. CUVSSEMENT OE L'OBJET OE LA DEMANDE 

CIB 7 G06N5/04 



Seton la classification intemationate des brevets (CIB) ou a la fois selon la dasstfication nationale et la CIB 
B. DOMAINES SUR LESQUELS LA RECHERCHE A PORTE 



Documentation minimaie consultee (systeme de classification suivi des symbdes de dassement) 

CIB 7 G06F G06N GOIV 



Documentation consultee autre que ta documentation minimaie dans ta mesure ou ces documents retevent des domaines sur lesquels a port6 la recherche 



Base de dormees electronlque corisultee au cours de ta recherche Internationale (nom de la base de donh^. et si realisable, termes de recherche utilises) 

EPO-Internal 



C. DOCUMEMTS CONSJDERES COMME PERTINENTS 



Categorie Identification des documents cites, avec, ie cas echeant, Tindication des passages pertinents 



no. des revendtcations visees 



CHU CH -K P ET AL: "FIRST BREAK 
REFRACTION EVENT PICKING USING FUZZY LOGIC 
SYSTEMS" 

IEEE TRANSACTIONS ON FUZZY SYSTEMS, US, IEEE 
INC. NEW YORK, 
vol . 2, no. 4, 

1 novembre 1994 (1994-11-01), pages 
255-266, XP000469682 
ISSN: 1063-6706 
abrege 

page 257, colonne de gauche, ligne 45 
-page 258, colonne de drolte, ligne 36 
page 265, colonne de gauche, ligne 15 - 
ligne 37 

EP 0 468 229 A (HNC INC) 

29 Janvier 1992 (1992-01-29) 

colonne 13, ligne 15 -colonne 15, ligne 9 



1,4,6 



1.4 



Voiria 



suite du cadre C pour la fin de la tiste des documents 



|){ j Les documents de families de brevets sont tndiqu^s en anr>exe 



* Categories spectales de documents dt^: 

*A' document definissant T^tat general de la technique, non 
considere comme particulierement pertinent 

'E' document anterieur, mais publie a la date de depot international 
ou apres cette date 

"L' document pouvant jeter un doute sur une revendlcation de 
priorite ou cite pour determiner la date de put>lication d*une 
autre citation ou pour une raison speciale (telle qu'indiquee) 

'0" document se referant a une divulgation orale, a un usage, a 
une exposition ou tous autres moyens 

"P* document publie avant la date de depot international, mats 
post^eurement d la date de priority revendlquee 



T* document uJtdrieur pubDe apres la date de depot international ou la 
date de priorite et n'appjartenenant pas a I'etat de la 
technique pertinent, mais cite pour comprendre le principe 
ou la th^e oonstituant la base de rinvention 

*X' document particufi^rement pertir>ent: Tinven tion revendlquee ne peut 
dtre consideree comme nouvelle ou comme impliquant une activite 
inventive par rapport au document considere isoi^ent 

*Y' document particulierement pertinent: I'inven tion revendlquee 
ne peut dtre consideree comme impliquant une activite inventive 
lorsque (e document est associe a un ou plusieuis autres 
documents dcf meme nature, cette combinatson etant ^vidente 
pour une peisorme du metier 

'&* document qui fait partle de la m^e famille de brevets 



Date a laquelle la recherche Internationale a ete effectivement achevee 



5 septembre 2000 



Date dexpedrtion du present rapport de recherche Internationale 



12/09/2000 



Nom et adresse postale de I'administration chargee de la recherche intemationaie 

Office Europeen des Brevets. P.B. 5818 Patentlaan 2 

NL-2280HV Rijswijk 

Tel. (+31-70) 340-2040. Tx. 31 651 epo nl. 

Fax: (♦31-70) 340-0016 



Fonctionnaire autorise 



Kingma, Y 



Forrmiairs PCTrtSA/210 (douxieme feuille) OuiHet 1992) 



RAPPORT DE "HERCHE INTERNATIONALE 



.6nit ' IntoniAtkMUilo No 

1 PCT/FR 00/01804 


C.(5ultB) DCX:UMENTS CONSIDERES COMME PERTtNEKTS 


Cat^rle ' 


tdenttrication dM documents clt6s, avecle cas M6ant, iMndleatlondes passages pertinents 


no. des revendealions vistes 


A 


MAGOTRA N ET AL: "SEISMIC SIGNAL 

DISCRIMINATION USING ADAPTIVE SYSTEM 

PARAMETERS" , PROCEEDINGS OF THE MIDWEST 

SYMPOSIUM ON CIRCUITS AND SYSTEMS, US, NEW 

YORK, IEEE, VOL. CONF. 33, PAGE(S) 84-87 

XP000280121 

1e document en entier 


1,4-6 



Foimulaife PCT/ISA«10 (siits do la douxjAfrw teuiB«) QuiUet 1992) 



RAPPORT DE RECHERC INTERNATIONALE 

Renaelgnements relatlfs au* mombresde families de brevots 



Denii b.».natlonato No 

PCT/FR 00/01804 



Document brevet cite 
au rapport de recherche 



Datede 
publication 



Membre{s) de la 
famille de brevet(s) 



Date de 
publication 



EP 0468229 



29-01-1992 



JP 
US 



5151188 A 
5398300 A 



18-06-1993 
14-03-1995 



RxTTulaire PCT/ISAr'210 (armaxe tamilesde brftv«ts) (juiliet 1902} 



10/0093^5 

B1 

Projet (NON destln6 ^ la presentation) imprim6 le 17.07.2000 10:28:37 AM 



0 

0-1 


Reserve a Toffice r^cepteur 
Demande Internationale No. 




0-2 


Date du depot international 




0-3 


Nom de Toffice recepteur et "Demande 
internationale PCT 








0-4 

0-4-1 


Formulaire-PCT/RO/101 Requite 
PCT 

Prepare avec 


PCT-EASY Version 2 . 90 
(mis a jour 10.05.2000) 


0-5 


Petition 

Le soussigne requiert que ia presente 
demande internationale soit trait^e 
conform6ment au Trait6 de cooperation 
en matiere de brevets 




0-6 


Office recepteur (choisi par le 
deposant) 


Ins ti tut national de la propriete ./ 
industrielle (France) (RO/FR) 


0-7 


Reference du dossier du deposant ou 
du mandataire 


B 13340.3 DB / 


1 


Titre de I'invention 


SYSTEME D * INTELLIGENCE ARTIFICXELI*E Jr'UUK 
LA CLASSIFICATION D » EVENEMENTS , D * OB JETS 
OU DE SITUATIONS A PARTIR DE SIGNAUX ET 
DE PARAMETRES DISCRIMINANTS ISSUS DE 
MODELES 


II 

11-1 
11-2 
IM 
11-5 

11-6 
11-7 
11-8 

II-9 


Deposant 

Cette person ne est : 
Deposant pour 
Nom 

Case* 

Nationalite (nom de I'Etat) 
Residence (nom de I'Etat) 
No. de telephone 

No de telecopieur: 


Deposant seulement 

Tous les Etats designes sauf US 

COMMISSARIAT A L ' ENERGIE ATOMIQUE 

rue de la Federation 
F-75752 PARIS 15eme 
France 
FR 
FR 

01 69 08 82 93 
01 69 08 82 92 


111-1 

111 1.1 
III- 1 - 1 

IIM-2 

JIM -4 

III-1-5 

III-1-6 
III-1-7 


Deposant et/ou inventeur 
Cette personne est : 

Deposant pour 

Nom (NOM DE FAMILLE. Pr6nom) 
Adresse: 

Nationalite (nom de I'Etat) 
Residence (nom de I'Etat) 


Deposant et inventeur 

US seulement 

MULLER, Jean-Denis 

5 allee de la Garenne 

F- 7 8120 CLAIREFONTAINE-EN-YVELYNES 

France 

FR 

FR 
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111-2 

m-2-1 


Deposant et/ou invent ur 
Cette personne est : 


Deposant et inventeur 




111-2-2 


Deposant pour 


US seulement 




II 1-2-4 


Norn (NOM DE FAMILLE. Pr6nom) 


MULLER CARCELES , Stephanie 




III-2-5 


Adresse: 


5 allee de la Garenne 

F-78120 CIAIREFONTAINE-EN-YVELYNES 
France 




III-2-6 


Nationalite (nom de I'Etat) 


FR 




III-2-7 


Residence (nom de I'Etat) 


FR 




IV-I 

IV-1-1 


Mandataire ; Reprdsentant commun 
ou adresse pour la correspondance. 

La personne nomm6e ci-dessous est/a 
ete designee pour agir au nom du ou 
des deposants aupres des autorites 
Internationales competentes, comme 
Nom (NOM DE FAMILLE. Prenom) 


mandataire 
LEHU , Jean 




IV-1-2 


Adresse: 


c/o BREVATOME 

3 / rue du Docteur Lancer eaux 
F~75008 PARIS > 
France 




IV-1-3 


No. de telephone 


01 53 83 94 00 




IV-1-4 


No de telecopieur: 


01 45 63 83 33 




IV-I-S 


Courrler ^lectroniquie: 


spibrev@ easynet . f r 




V 


Designation d'Etats 




V-1 


Brevet regional 

(d'autres formes de protection ou de 
traitement, le cas echeant, sont 
specifiees entre parentheses pour les 
Etats designes concernes) 


EP: AT BE CH&LI CY DE DK ES FI FR GB GR 
IE IT LU MC Nli PT SE et tout autre Etat 
qui est un Etat contractant de la 
Convention sur le brevet europeen et du 






PCT 




V.2 


Brevet national 

(d'autres formes de protection ou de 
traitement, le cas echeant, sont 
specifiees entre parentheses pour les 
Etats designes concernes) 


CA JP US 


V-5 


Declaration concernant les 
designations de precaution 
Outre les designations faites sous les 
rubriques V-1 , V-2 et V-3, le deposant 
fait aussi, conform^ment a la regie 
4.9.b), toutes les designations qui 
seraient autorisees en vertu du PCT, e 

indiquee(s) dans la rubrique V-6 
ci-dessous.Le deposant declare que ces 
designations additionnelles sont faites 
sous reserve de confirmation et que 
toute designation qui n'est pas 
confirmee avant I'expiration d'un delai de 
15 mois a compter de la date de priorite 
sera consideree comme retiree par le 
deposant d I'expiration de ce deiai. 




V-6 


Exclusion(s) des designations de 
precaution 


NEANT 
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VI-1 
VI-1-1 


Revendication de priorite d'une 
demande Rationale anterieure 


•^^f'i 1 1 A't- 1 QQQ 




.07.1999) 


VI-1 -2 












VI-1 -3 


Pays 


31 X\ 








VIM 


Administration chargee de ta 
recherche Internationale choisie 


Office europeen 
(ISA/EP) 


des 


brevets (OEB) 


VI 1-2 


Demande d'utilisation des resultats 
d'une recherche ant^rieure; mention 
de cette recherche 

Date 


26 mai 2000 (26. 


05. 


2000) 




iNumcro 


FA 581821 








V/l 1-9 -"^ 




EP 








VIII 


Bordereau 


Nombre de feuilles 


Do5sier(s) electronique(s) joint(s) 


V Ml- 1 




4 




VI!l-2 


Description 


35 




VIII-3 


Revendications 


2 




VI 11-4 


Abrege 


1 


abr-13340db . txt 


vni-5 


Dessins 


12 




VIII-7 


TOTAL 


54 - 




Elements joints 


Document(s) papier joint(s) 


Dossier(s) electronique(s) joint(s) 


VI!l-8 


Feuitle de calcul des taxes 






VIII-10 


Copie du pouvoir general 


reference n** PG 
07085 








VIII-16 


Disquette PCT-EASY 




dlsque-tte 


\/lll-17 

VIII- 1 / 


nunc y^icwidci/ . 


liis^be des 
manda t a i r e s 
BREVATOME. 








Vill-18 


Figure des dessins qui doit 
accompagner I'abrege 


4 


VIII-19 


Langue de depot de la demande 
Internationale 


f ranqiais 


IX-I 


Signature du deposant ou du 
mandataire 










IX-1-1 


Norn (NOM DE FAMILLE. Prenom) 


LEHU , Jean ^ 









RESERVE A L'OFFICE RECEPTEUR 



10-1 


Date effective de reception des 
pieces supposees constituer la 
demande Internationale 




10-2 


Dessins: 




10-2-1 


Re^us 




10-2-2 


non re^us " 
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10-3 


Date effective de reception, rectifiee 
en raison de la reception ulterieure, 
mats dans les delais, de documents 
ou de dessins completant ce qui est 
suppose constituer la demande 
Internationale: 




10-4 


Date de reception, dans les delats, 
des corrections demanddes selon 
Particle 11.2) du PCT 




10-5 


Administration chargee de la 
recherche Internationale 


ISA/EP 


10-6 


Transmission de la cople de 
recherche differee jusqu'au paiement 
de ia taxe de recherche 






RESERVE AU BUREAU INTERNATIONAL 


11-1 


Date de reception de I'exemplaire 
original par le Bureau international 
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(Cette feuille ne fait pas partie de la demande intemationale ni ne compte comme une feuitle de celle-d) 



0 

0-1 


Reserve a I'office recepteur 
Demande intemationale No. 




0-2 


Timbre a date de ('office recepteur 






0-4 
0-4-1 


Formulaire - PCT/RO/101 (Annexe) 
Feuille de calcul des taxes PCT 
Prepare avec 


PCT-EASY Version 2 . 90 
(mis a jour 10.05.2000) 


0-9 


Reference du dossier du deposant ou 
du mandataire 


B 13340.3 DB 


2 


Deposant 


COMMISSARIAT A L * ENERGIE AT 
al. 


OMIQUE, et 


12 

12-1 
12-2 
12-3 

12^ 
12-5 
12-6 
12-7 
12-8 

12-9 

12-10 

12-11 

12-12 
. 12-13 


Calcul des taxes prescrites 


Montant total des 
taxes/muttiplicateur 


Montant total (FRF) 




Taxe de transmission T 




400 




Taxe de recherche S 




6 198.79 




Taxe intemationale 
Taxe de base 
(30 premiferes feuilles) b1 


2 682.86 




Feuilles suivantes 


24 


Montant additionnel (X) 


59.04 


Montant total additionnel b2 


1 416.96 


b1 + b2 = B 


4 099.82 


Taxes de designation 
Nombre de designations indiquees 
dans la demande intemationale 


4 


Number of designation fees 
payable (maximum 8) 


4 


Montant de la taxe de (X) 
designation 


577.24 


Montant total des taxes de D 

designation 


2 308.96 


Reduction de taxe POT-EASY . R 


-826.51 


Montant total de la taxe , 1 
intemationale (B+D-R) 


CO 


5 582.27 




12-17 


TOTAL DES TAXES DUES (T+S+l+P) 




12 181.06 




12-19 


Mode de paiement 


autorisation de debiter un compte de 
depot 


12-20 


Instructions concernant le compte de 
depot 

Loffice recepteur: 


Institut national de la propriete 
industrielle (France) (RO/FR) 


12-20-1 


est autorise a debiter mon compte de 
depot du total des taxes indique 
ci-dessus 
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12-20-2 


est autorise d debiter mon compte de 
depot de tout montant manquant. ou k le 
crediter de tout excedent, dans le 
paiement du total des taxes indique 
ct-dessus 




12-21 


Compte de d^pdt No. 


024 


12-22 


Date 


28 juin 2000 (28.06.2000) 


12-23 


Nom et signature 


DUPONT Catherine,.-..--^ 




MESSAGES DE VALIDATION ET REMARQUES 


13-2-1 


Messages de validation 
Requete ' 


Vert? 

Le titre de 1 ' Invention doit etxe bref 
et precis . Priere de verifier. 


13-2-2 


Messages de validation 
Etats 


Vert? 

11 est possible d'effectner davantage de 



pas ete designes : AP:( GH, GM, KE, LS , 





MZ, 


SD, 


SL, 


SZ, 


TZ, 


UG, 


ZW) 


; EA: ( 


AM, 


AZ, 


BY,. 


KG, 


KZ, 


MD, 


RU, 


TJ, 


TM) ; 


OA: 


( BF 


, BJ 


, CF 


, CG 


/ci 


, CM 


, GA 


, GN, GW, 


ML, 


MR, 


NE, 


SN, 


TD, 


TG) 


; A£ 


, AG 


, AL, AM, 


AT, 


AU, 


AZ, 


BA, 


BB, 


BG, 


BR, 


BY, 


BZ, CH, 


LI, 


CN, 


CR, 


CU, 


CZ, 


DE, 


DK, 


DM, 


DZ, EE, 


ES, 


FI, 


GB, 


GD, 


GE, 


GH, 


GSA, 


HR, 


HU, ID, 


IL, 


IN, 


IS, 


KE, 


KG, 


KP, 


KR, 


KZ, 


LC, LK, 


LR, 


LS, 


LT, 


LU, 


LV, 


MA, 


MD, 


MG, 


MK, MN, 


MW, 


MX, 


MZ, 


NO, 


NZ, 


PL, 


PT, 


RO, 


RU , SD , 


SE, 


SG, 


SI, 


SK, 


SL, 


TJ, 


TM, 


TR, 


TT , TZ , 


UA, 


UG, 


uz. 


VN, 


YU, 


ZA, 


ZW. 


Priere de 



verifier . 



Messages de validation 
Bordereau 



Messages de validation 
Paiement 



Vert? 

Priorite 1 : Le document de priori te 
n'est pas joint. (Le deposant doit le 
fournir dans un delai de 16 mois a 
compter de la date de priorite la plus 
ancienne r e vendiq[uee ) '_ 



Vert? 

Priere de verifier que vous avez bien un 
compte de depot aupres de 1 • office 
recepteur , • 



